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ABSTRACT 


* 

The  following  paper  studies  the  possibilities  of  using  high-speed  flywheels 
with  magnetic  bearings  in  communications  spacecraft.  Hardware  point  designs 
have  been  completed  for  a very  low  drag  bearing  and  a very  efficient,  compatible 
motor  generator.  Using  the  results  of  these  designs,  a tradeoff  study  compares 

if 

these  flywheels  with  combinations  of  NiCd  batteries,  NiH2  batteries,  and  con- 
ventional low-speed  momentum  wheels. 

Section  I presents  the  tradeoff  studies  for  flywheels  used  for  energy 

storage  only,  for  attitude  control  only,  and  for  combined  energy  storage  and 

attitude  control.  The  results  show  that  such  wheels  look  excellent  for 

attitude  control  alone  (as  high-speed,  high-reliability  momentum  wheels) 

and  attractive  for  combined  power  and  attitude  control,  but  not  much  better 

than  NiH1  for  power  alone,  /f  — 

^ v 

Section  II  reviews  the  technology  of  flywheel  components,  showing  that 
rotors  are  the  biggest  unsolved  problem  but  can  be  built  now  to  lower  energy 
densities. 

Section  III  analyzes  two  types  of  anisotropic  rotors  to  allow  design 
calculations  for  a given  set  of  material  properties.  These  types  are  the 
orthotropic  composite  rotor  and  the  bare  filament  rotor. 

Section  IV  documents  the  design,  development,  and  testing  of  the  Lincoln 
Laboratory  magnetic  bearing.  This  bearing  can  be  used  in  momentum  wheels 
which  are  not  electronically  gimballed,  or  in  energy  storage  wheels. 

Section  V documents  the  development  of  the  Lincoln  Laboratory  motor 
generator,  which  allows  very  high  efficiency  at  high  speeds. 

Section  VI  contains  overall  conclusions  of  this  study  concerning  space- 
craft applications  of  high-speed  flywheels. 
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I.  SYSTEM  ASPECTS  OF  FLYWHEELS 
A.  Flywheel  Control 

A flywheel  system  will  be  considered  in  the  context  of  a synchronous 
comraunicat ion  satellite  mission.  Assume  a 1 KW  payload  and  an  attitude 
control  system  with  required  stability  of  ± .01°  with  dominant  solar  torques 
on  the  order  of  50  microfoot  pounds  sinusoidal  component. 

The  Lincoln  Laboratory  (LL)  flywheel  unit  concept  briefly  consists  of  a 
rotor,  shaft,  magnetic  bearings,  touchdown  bearings,  motor-generator,  vacuum 
housing,  and  miscellaneous  small  parts.  To  ease  unit  integration  for  a first 
cut  we  have  selected  a small  shaft  with  essentially  zero  internal  rotor  diam- 
eter, although  this  costs  a factor  of  two  in  rotor  weight  compared  to  the  limit 
of  ID/OD  -*■  1.  We  will  assume  two  levels  of  rotor  technology:  a 34  watt-hour 
per  pound,  40  K RPM  rotor,  probably  achievable  with  present  technology  and 
minor  improvements  in  fabrication  techniques;  and  a 54  watt-hour  per  pound 
50  K RPM  rotor  of  the  same  size  (20-in.  dia.),  probably  the  ultimate  for 
present  materials  and  the  given  geometry.  Vacuum  containment  was  chosen 
rather  than  lightweight  covers  for  flight  to  provide  complete  testing  of  a 
fully  assembled  unit  and  to  add  bearing  mount  rigidity  and  small  fragment 
conf inement . 

These  units,  denoted  LL-30  and  LL-50,  are  estimated  to  weigh  respectively 
45  and  50  pounds.  They  will  be  considered  first  as  purely  energy  storage 
devices  in  competition  with  batteries,  then  as  momentum  wheels  in  competition 
with  presently  available  and  current  prototype  designs,  and  finally  as  com- 
bined energy  and  angular  momentum  storage  devices. 

The  essential  relationships  between  wheel  speed,  energy,  momentum,  and 
mass  are  given  below  and  graphed  in  Fig.  1-1.  A disc  rotor  with  inner  radius 
Rj  and  outer  radius  R^,  mass  density  p and  thickness  h has 


mass 

M = 

nph 

2 2 
(Ri  - V 

(1-1) 

moment  of  inertia 

I = 

1/2 

M (R2  + R2) 

I o 

(1-2) 

angular  momentum 

H = 

IU) 

(1-3) 

j_18  8 190S41 


H'  « TOTAL  NORMALIZED 
MOMENTUM 

E'  * TOTAL  NORMALIZED 
ENERGY  (2  wheels) 

u,,u2  --  WHEEL  SPEEDS 
uM  = MAXIMUM  SPEED 


SHADED  AREA  IS  ACCESSABLE 
TO  COMBINED  ENERGY  AND 
ATTITUDE  CONTROL  SYSTEM 


Fig.  1-1.  Flywheel  parameter  relationship. 


energy  E = 1/2  lo^  = 1/2  Ho;  (1-4) 

energy/mass  E/M  = 1/4  (R^  + R~)  o_"  (1-5) 

2 "> 

momentum /mass  H/M  = 1/2  (R  + R“)  u (1-6) 

I o 

For  our  case  RT  = o and  for  normalized  variables  Ui'  = to/co  where  Oj  = maximum 
I ^ m m 

speed.  Total  energy  of  two  wheels  E'  = E/Iti/,  H'  - H/lrn^. 

E'  = 1/2  (ui'2  + u'2)  (1-7) 

H'  = - w'  (1-8) 


A big  concern  with  any  energy  wheel  system  on  a satellite  will  be  the 
large  angular  momentum  (H)  stored  in  their  energy  system,  and  its  potential 
effects  on  attitude  control.  A basic  ground  rule  of  any  such  system  must  be 
that  no  likely  single  point  failure  shall  cause  total  loss  of  the  satellite. 
This  implies  that  the  following  failure  modes  must  not  be  very  likely. 

1.  Rotor  breaking  catastrophically,  throwing  high-energy  fragments. 

2.  Suspension  totally  lost  at  high  speed,  causing  rotor  failure. 

3.  Wheel  overspeed. 

On  a less  stringent  level,  to  avoid  a recoverable  tumbling  of  the  satel- 
lite, the  following  would  be  desirable: 

1.  ACS  automatic  rate-dumping  mode  not  dependent  on  any  one  wheel. 

2.  Wheel  speed  control  redundancy. 

3.  Simple,  reliable  suspension  with  redundancy  and  recovery 
capability. 

4.  A touchdown  bearing  system  capable  of  at  least  100  touchdowns 
and  100  hour  continuous  running  at  design  speed,  perhaps  with 
increased  power  requirements. 

5.  Simple  energy  conversion  from  generator  to  suspension,  includ- 
ing touchdown  recovery. 
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The  result  of  these  conditions  would  then  be  a flywheel  energy  system 
which  would  not  be  likely  to  disturb  the  attitude  control  of  the  satellite 
in  case  of  component  failure,  and  would  not  result  in  total  loss  of  satellite 
even  if  it  did  so. 

The  implications  for  construction  of  the  flywheel  are  as  follows: 

1.  High  confidence  in  the  rotor  not  to  break  catastrophically. 

2.  Redundant  overspeed  control. 

3.  A radially  passive  suspension  system,  the  simplest  and  most 
reliable  choice,  is  very  attractive  unless  active  axial 
suspension  can  be  made  fully  redundant  without  compromising 
simplicity  and  reliability. 

4.  Redundant  axial  suspension,  possibly  with  ground  intervention 
by  command . 

5.  A reliable  high-speed  touchdown  bearing  system. 

6.  Redundant  motor-generator  windings  if  the  wheel  energy 
storage  function  must  be  recovered. 

7.  Efficient,  low-power  suspension. 

In  looking  at  the  actual  realization  of  such  a wheel  system,  several 
points  became  clear: 

1.  Most  realizations  of  active  suspension  electronics  include 
several  failure  modes  which  result  in  touchdown  at  high  speed 
with  reduced  chance  of  automatic  recovery,  but  with  ground- 
commanded  recovery  capability.  Very  careful  design  in  this 
area  may  result  in  a low  probability  of  touchdown.. 

2.  A radially  passive,  axially  active  (RPAA)  magnetic  suspension 
with  a V-groove  mechanical  touchdown  skid  on  the  ball  bearing 
races  gives  a simple  caging  mechanism  for  launch;  simply 
remove  bearing  power  or  bottom  the  wheel  actively. 

3.  The  RPAA  magnetic  suspension  results  in  a penalty  ir  spinning 
energy  losses  due  to  eddy  current  drag.  There  is  a great 
deal  of  work  to  be  dor.e  in  the  area  to  minimize  the  drag. 

Wheel  balancing  is  important  to  minimize  this  loss.  Low 
stiffness  ratio  (see  below)  is  also  important. 
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4.  The  RPAA  suspension  requires  a broadband  audio  frequency  servo- 

loop  design  with  possible  complications  due  to  acoustic  resonances 

of  the  structure.  The  loop  bandwidth  is  approximately  equal  to 

the  time  constant  A'. , /M  where  M is  the  rotor  mass  and  K.  is  the 

A A 

axial  negative  spring  constant  of  the  passive  magnetic  suspension. 

This  loop  must  be  designed  into  the  structure  of  the  satellite, 

and  structural  modes  included  in  the  loop  design. 

5.  The  axial  spring  constant  is  about  five  to  ten  times  the 

positive  radial  spring  constant  K which  passively  centers 

K 

the  rotor.  For  1G  operation  and  the  required  cross-axis  rates 
cuc , this  stiffness  must  be  large  enough  to  avoid  touchdown. 

The  minimization  of  stiffness  ratio  is  a development  area  with 
big  payoff. 

6.  The  touchdown  system  is  most  reliable  with  small  (several  mils 
to  ten  mils)  clearances.  This  implies  clearance  S less  than 
some  maximum,  limiting  motions  for  a fixed  disturbing  force 

or  torque,  and  so  requiring  a high  stiffness  from  the  magnetic 
bearings  and  a wide-band  servo  problem  with  acoustics  to 


complicate  matters. 

2 

Cross  axis  torque: 

(kr  l"ka  -l>  l > “c  H 

(1-9) 

1G  support: 

KD  S > MG 

i\ 

(1-10) 

Stiffness  Ratio: 

K > N Kn  5 < N < 10 

A R 

(1-11) 

Plugging  in,  floQp 

V 

ro|»— 

=«  1 

CO  | O 1 

(1-12) 

For  S = 10  mils,  N = 8 this  gives  f,  > 89  Hz,  a range  in  which  careful 

Loop  ° 

design  is  required  to  avoid  acoustics.  Also,  sensing  elements  at  this  band- 
width must  be  carefully  designed  to  avoid  noise  pickup  and  phase  shift. 

Notice  that  the  cross-axis  torque  equation  (1)  puts  a limit,  not  only  on 
minimum  K^,  but  on  the  aspect  ratio  of  the  fringing  rings  D/L,  where  D = 
diameter  of  ring  and  L = spacing  along  the  axis  between  bearings.  For  good 
cross-axis  torque  we  need  D/L  <<  K^/K^.  This  determines  the  system  axial 
length . 
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7.  Nutation  damping  may  be  a serious  limitation  to  wheel  speed. 

For  damping  at  nutation  frequencies  in  the  suspension  system, 
wide  loop  bandwidths  (several  to  10  times  rotational  rate) 
are  needed  and  wheel  inbalance  couples  into  the  suspension  to 
produce  losses.  This  area  needs  work  for  an  imbalance-tolerant 
system;  composite  wheels  tend  to  suffer  from  balance  problems. 

With  these  matters  in  mind,  the  LL  flywheel  unit  has  been  designed. 

The  following  sections  use  it  as  a basis  for  system  trade-off  studies. 

B . Energy  Storage  Only 

The  basic  energy  storage  unit  for  a flywheel  system  is  a pair  of  wheels 
controlled  to  give  zero  total  angular  momentum.  For  minimum  u)'  = 0.4,  Eq . (1-7) 
gives  an  accessible  energy  AE'  = 0.84.  This  limit  occurs  because  the  motor- 
generator  is  effectively  torque-limited.  At  constant  H,  the  torque  required 
for  a power  level  P is 

T ° 1 It  = P/(wl  + 'V  (1-13) 

This  is  maximum  at  minimum  oj. 

The  characteristics  of  the  flywheel  unit  are  listed  in  Table  1-1.  Note 
that,  for  power  only,  the  flywheel  accessible  energy  to  weight  ratio  for  the 
rotating  unit  is  a bit  more  than  a third  that  of  the  rotor  alone.  Then,  to 
convert  the  power  out  to  usable  DC,  the  electrical  efficiencies  and  weights 
of  machine,  electronics,  and  other  components  are  estimated  in  Table  1-2.  This 
costs  another  factor  of  two  in  available  energy  to  system  weight.  However, 
the  final  figures  are  quite  respectable. 

The  competitor  battery  systems  fall  into  several  classes.  Batteries  may 
be  NiCd  or  NiH2«  They  may  be  operated  at  various  depth  of  discharge  (DOD) 
levels,  with  statistics  for  various  resulting  lifetimes.  They  may  be  operated 
through  charger  and  discharger  electronics  for  a fully  regulated  bus,  but  this 
eliminates  the  inherent  advantage  of  a battery  as  a DC  voltage  source.  The 
tradeoff  study  shown  in  Table  1-3  presents  the  options.  Option  I represents 
a NiCd  lightweight  battery  with  50%  battery-level  redundancy,  60%  DOD,  with 
charger  but  no  discharge  converter,  so  that  the  bus  floats  at  the  battery  out- 
put voltage  level  during  eclipse.  Option  II  is  a NiH^  battery  on  a fully 


6 


TABLE  1-1 


SINGLE  FLYWHEEL  CHARACTERISTICS 


Rotor  Type  LL-30 

Rotor,  10-in.  radius,  lbs.  (Kg)  20  (9.08) 

Motor-generator,  lbs.  5 

Bearings,  lbs.  10 

V iutn  housing,  lbs.  5 

Shaft,  lbs.  3 

Miscellaneous  lbs.  2 

Total  weight,  lbs.  45 

Speed  range,  RPM  16K-40K 

Max.  Rotor  E/M,  w-h/lb.  34 

Max.  Stored  E,  w-h  (MJ)  680  (2.45) 

Available  E,  w-h,  H'  *=  0.1  476 

Available  E,  w-h,  H " *=  0 571 

Available  E/M,  H ' - 0 12.7 

1.2  hr.  Discharge  Rate,  H'  * 0 475 

2 

Moment  of  Inertia,  Kg-m  0.29 

H,  Angular  Momentum,  ft. -lb. -sec.  (N-M-S)  867  (1188) 


LL-50 

20  (9.08) 
10 
10 
5 
3 
2 

50 

20K-50K 

54 

1080  (5.89) 
756 
907 
18.1 
755 
0.29 

1084  (1485) 
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TABLE  1- 


b 


2-WHEEL  POWER-ONLY  SYSTEM 
DISCHARGING  CHARACTERISTICS  IN  ECLIPSE 


Output  powe:  from  wheels,  W 
Generating  machine  efficiency  % 

Generating  bridge  efficiency  % 

Drag  power  from  bus,  W 

Suspension  and  control  power  from  bus,  W 

Output  power  to  load 

Total  generating  efficiency  % 

Bridge  weight,  lbs.,  .035  lb. /watt 
Control  electric  weight,  lb.,  .005  lb. /watt 
Suspension  electric  weight,  lb. 

Total  electric  weight,  lb. 

Total  system  weight 

Total  system  output  energy  (w-h)  1.2  hr. 
Total  system  E/M,  w-h/lb. 


LL-30  LL-50 


950 

96 

92.5 
20 

5 

819 

86 

76 

3.6 

5.0 

34.6 
124.6 
983 

7.9 


1510 

96 

92.5 

20 

5 

1316 

87 

42 

6.0 

5.0 

53.0 

153.0 

1579 

10.3 


20 

5 

182 

94.5 

89.5 
240 

76 
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CHARGING  CHARACTERISTICS  IN  SUNLIGHT,  10  HOUR  RATE 


Average  drag  losses,  W 20 

Suspension  and  control  power,  W 5 

Wheel  input  power,  W 114 

Motoring  efficiency,  motor,  7 94.5 

Motoring  efficiency,  bridge,  % 89.5 

Total  input  power,  W 160 

Charging  efficiency,  % 71 


TABLE  1-3 


BATTERY  SYSTEM  COMPARISONS 
1.0  Kw  System 


Option  1:  NiCd  batteries,  50%  battery  redundant,  with  charger  but  no 

discharges,  only  diode  isolation  and  protection,  6.8  w-h/lb., 
60%  DOD. 


Option  1 1 : NiH^  batteries,  20%  cell  redundant,  with  charges  and  dis- 

charger for  fully  regulated  bus,  75%  DOD,  14  w-h/lb. 

Option  III:  NiH?  batteries,  20%  cell  redundant,  with  charger  but  only 

output  isolation  and  protection,  75%  DOD,  17  w-h/lb. 


I 

II 

III 

Control  power,  W 

5 

5 

5 

Discharge  efficiency 

.95 

.85 

.95 

Total  discharge  power,  W 

1057 

1181 

1057 

Stored  E,  w-h 

1263 

1417 

1263 

Charger  efficiency,  c/10 

.8 

.8 

.8 

Heater  power,  W 

40 

40 

40 

Charging  input  power,  W 

197 

217 

197 

Charger  and  heater  weight,  lb. 

(.11  lb./W) 

21.7 

23.9 

21.7 

Discharger  weight,  lb. 

(.05  or  .025  lb/W) 

26 

59 

26 

Total  electronics  weight,  lb. 

48 

83 

48 

Battery  weight  with  redundancy,  lb. 

262 

101 

90 
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regulated  bus,  with  20%  cell  level  redundancy  and  75%  DOD.  Option  III,  the 
lightest,  is  Option  II  without  the  discharge  converter. 

The  resulting  figures  in  Table  1-4  show  system  weight  at  1 KW  load, 
system  E/M,  input  charging  load  on  the  solar  arrays,  extra  resulting  solar 
array  weight,  and  total  weight. 

Based  on  these  figures,  it  is  easy  to  see  that  either  Nil^  or  flywheel 
storage  is  much  lighter  than  NiCd  batteries.  Both  are  developmental,  but 
NiH^  is  far  more  advanced  and  no  heavier.  Because  of  this,  there  does  not 
appear  to  be  good  reason  to  pursue  flywheels  purely  for  energy  storage. 

Even  the  LL-50  system  with  its  higher  energy  density  does  not  look  much 
better  than  nickel-hydrogen  batteries. 

Note  also  that  the  conditioning  electronics  become  significant  sources 
of  power  loss  and  system  weight  in  both  Option  II  and  LL-50.  This  implies 
that  the  inherent  DC  voltage  source  nature  of  the  chemical  battery  seems  to 
give  lightweight  Nil^  or  similar  systems  an  ultimate  edge  over  kinetic  storage. 

C . Attitude  Control  Only 

A high  energy  density  flywheel  should  make  a very  lightweight  momentum 
wheel.  At  first  glance,  the  optimization  would  favor  a larger  diameter  and 
lower  speed  for  best  momentum  density.  However,  this  tradeoff  very  quickly 
is  limited  by  the  size  of  the  wheel  and  the  form  factor  of  the  enclosure.  In 
fact,  simply  by  reducing  the  speed  of  an  energy  wheel  by  a healthy  factor  to 
reduce  the  stresses  gives  a rather  attractive  momentum  storage  device. 

How  much  angular  momentum  is  enough?  Several  criteria  could  be  used. 

For  a LES-8/9  type  biased  momentum  system,  the  maximum  solar  torque  acting 

over  the  6 hours  of  uncontrolled  yaw  should  produce  a negligible  yaw  rotation. 

Taking  worst  case  numbers  of  50  microfoot-pounds  and  yaw  e = 1.2°  gives  a bias 

requirement  of  H - 2T/uo£  “ 68  f t . -lb. -secs . , a moderate  figure,  for  daily 

u i ‘ 7.3  x 10  5 rad/sec. 
o 

Another  criterion  would  be  to  avoid  any  thrusting  due  to  the  daily  solar 
torque  variation.  For  a daily  U)Q  we  have 
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TABLE  1-4 


POWER- ONLY  WEIGHT  TRADEOFF 
1.2  KW-Hr.  Battery 


1 

n. 

III 

LL-30 

LL-50 

X 1.22 

X .760 

Total  system  weight,  lbs. 

310 

184 

138 

152 

116 

System  E/M,  W-hr./lbs.  usable 

3.87 

6.52 

8.69 

7.9 

10.3 

Input  charging  load,  W 

197 

217 

197 

195 

182 

Extra  solar  array  wt.,  .11  lb./W 

22 

24 

22 

21 

20 

Total  Weight  Penalty,  lbs. 

332 

208 

160 

173 

136 

11 


AH 


T Sin  Li)  t dt 
o o 


(1-14) 


For  Tq  = 50  pft.-lb.  this  gives  AH  » 1.37  ft. -lb. -sec.  If  this  is  to  tip  over 
the  wheel  by  no  more  than,  say,  1.2°,  this  gives  a momentum  requirement 

AH  -9 

— — £ 2 x 10  ^ radians  (1-15) 

or  H > 70  ft. -lb. -sec.  This  is  a moderately  large  angular  momentum,  as  deter- 
mined in  the  next  section.  To  avoid  tipping  over  by  0.12°  would  require  700 
f t . -lb . -sec . , certainly  a large  wheel,  about  the  size  of  the  energy  storage 
prototype  wheel  (45  lbs.)  using  advanced  flywheel  technology. 

One  criterion  of  angular  momentum  size  is  to  ask  whether,  if  all  the 

angular  momentum  of  the  wheel  were  suddenly  transferred  to  the  spacecraft,  the 

3 2 

craft  would  survive.  For  a spacecraft  moment  of  inertia  of  10  Kg-m  and  a 
maximum  recoverable  body  rate  of  1 RPM  (W  «*  .1  rad/sec)  this  would  allow 
H = 100  N-m-s  or  74  ft. -lb. -sec.  Another  criterion  would  be  whether  the 
angular  momentum  could  change  by,  say,  1%  and  not  seriously  interfere  with 

_3 

operations.  For  allowable  rates  of  .05°/sec.  (w  = 10  ) this  would  give 

H - 100  f t. -lb. -sec. , so  the  criteria  are  roughly  equivalent.  By  these 
criteria,  any  zero  bias  system  will  have  small  momentum,  and  the  minimum 
2.5  ft. -lb. -sec.  bias  is  also  small.  A biased  system  avoiding  daily  thrust 
is  marginal.  An  energy  storage  system  is  large. 

To  compare  attitude  control  capabilities  alone,  consider  four  possibil- 
ities. First,  consider  only  a 10  ft. -lb. -sec.  wheel  with  enough  extra  fuel  to 
last  10  years  thrusting  against  daily  solar  torque  variations.  Second,  consider 
four  10  ft. -lb. -sec.  reaction  wheels  with  no  bias  momentum.  Third,  consider  a 
100  ft. -lb. -sec.  wheel  with  conventional  construction.  Last,  consider  a LES-30 
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wheel  scaled  down  to  a 7-in.  radius  and  5-lb.  rotor.  The  system  weight  trade- 
off and  LES  weight  budget  is  shown  in  Table  1-5. 

To  arrive  at  fuel  penalty  weight,  we  have  used  a specific  impulse  of  gas 

of  100  secs.,  and  a moment  arm  of  2.5  ft.  Solar  torque  at  4T  / cu  for  T = 

' o o 

50  pft.-lb.  gives  3 ft. -lbs.  per  day  for  3650  days  totals  1.1  x 10  ft. -lb. -sec. 

0 4 

or  44  lbs.  of  gas.  Wheel  weights  are  taken  from  the  rule  M = 7H  which  gives 
a good  fit  to  Bendix  wheel  weights  over  a wide  range,  plotted  on  Fig.  1-2. 

Note  that  the  scaled-down  LL  flywheel  is  not  that  far  in  momentum  or 
weight  from  the  Teldix  MSR  momentum  ring,  which  was  an  iron  rotor  and  5-axis 
active  magnetic  bearings  with  a hollow  ring  geometry  to  achieve  high  momentum 
at  lower  speeds.  (See  Eq.  6:  if  = Rq,  then  w can  be  half  as  large  for  the 
same  H.) 

The  result  of  the  tradeoff  is  clear:  the  weight  and  complexity  of  the 
4-wheel  Option  II  is  poor.  The  weight  of  Option  I (small  wheel,  extra  gas)  Is 
slightly  higher  than  Option  III  (large  conventional  wheel)  but  both  are  far 
heavier  than  Option  IV,  the  advanced  flywheel.  The  message  is  that  lighter, 
longer-lasting  momentum  wheels  can  now  be  built  and  should  be  used. 

One  other  possibility  is  very  attractive  for  a high-speed  magnetically 
suspended  wheel.  Magnetic  gimballing  of  the  flywheel  rotor  is  possible  without 
adding  any  extra  hardware,  provided  that  the  offset  angle  is  small  enough.  The 
Teldix  wheel  with  all  active  suspension  can  do  this  simply  by  changing  the  off- 
set on  the  four  cross-axis  position  sensors.  A permanent  magnet  passive  sus- 
pension can  also  be  modified  to  do  this,  but  with  more  drag  in  the  offset 
(gimballed)  condition  perhaps.  Rotation  is  limited  by  bearing  gap  geometry 
to  0.5  to  1.2  degrees,  depending  on  the  design. 

To  allow  complete  3-axis  stabilization  (no  angle  at  all)  the  cross-axis 
momentum  must  equal  the  integrated  daily  solar  torque  over  a day.  From  our 
assumptions  and  Eq.  (10)  we  need  1.37  ft. -lb. -sec.  Allowing  1.2°  gimballing 
gives  a 70  ft. -lb. -sec.  wheel  (100  n-m-s)  which  can  be  achieved  with  the 
existing  Teldix  MDR  momentum  wheel.  An  improved  subscale  version  of  this 
now  has  fully  redundant  windings  for  all  bearings  and  motor,  providing  the 
needed  reliability  for  a long-duration  version.  Therefore,  the  MDR  type 


13 


TABLE  1-5 


GPSCS  ACS  TRADEOFF 


1 

Single  10  FPS 
Wheel 


Bias  momentum  - 10 

Ft . -lb . -sec . 

Weight,  1 wheel  20 

Weight,  all  wheels  20 

Weight  fuel  penalty  44 

Total  weight  50 

Structure  Penalty  Weight  0 

Elect,  weight  5 

Total  Weight  69 

Complexity  good 


15-lb.  100  FPS  Wheel  LL-30 

Rotor  4 OK  RPM,  R - 7”  M 
M/G 

Mag  bearing 
Housing 

Shaft  and  misc. 
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III 

IV 

4-10  FPS 

100  FPS 

LL-30 

Wheels 

Wheel 

Wheel 

0 

100 

100 

20 

50 

15 

80 

50 

15 

0 

0 

0 

80 

50 

15 

20 

0 

0 

20 

10 

10 

120 

60 

25 

poor 

good 

good 

Variety-Weight  Budget 

5.0  lb.  E/M  = 16  W_H/lb. 

2.0  lb. 

3.0  lb. 

3.0  lb. 

2.0  lb. 

15.0  lb. 


MASS  { lbs) 


1000 


100 


10 


4 


13  7 

|ia  \ i89i?]  * 


137 

“I 


40 


10 


I 


100  1000 
M (ft  - lb  jec ) 


Fig.  1-2. 


Momentum  wheel 


rotating  unit  tradeoffs. 
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MASS  Ikql 


wheel  at  70  ft. -lb. -sec.  with  redundant  windings  is  probably  the  best  candi- 
date for  a general-purpose  platform  with  3-axis  stabilization  and  low  fuel 
util iza t ion . 

The  unit  weighs  only  about  20  lbs.  and  has  no  wear-out  mechanism  which 
appears  worse  than  electronics  reliability  for  long-term  use  in  space.  The 
only  single-point  failure  would  be  a short  between  turns  in  the  motor  or 
bearing  coils,  so  these  should  be  carefully  inspected  and  tested  and  treated 
as  a mechanical  defect  rather  than  a random  electronics  failure.  The  unit 
would  be  able  to  operate  over  a much  wider  range  of  temperatures  than  previous 
ball-bearing  systems,  and  would  therefore  allow  a simplification  of  the 
satellite's  thermal  control. 

Table  1-6  gives  requirements  for  such  a wheel.  Figure  1-3  shows  construc- 
tion of  an  MDR  type  Teldix  wheel,  similar  to  that  proposed  here.  The  Teldix 
MDK  wheel  with  increased  gimballing  angle  and  full  redundancy  in  all  coils  and 
windings  represents  the  best  general-purpose  fully  stabilized  platform  ele- 
ment presently  near  complete  development. 

D . Flywheel  Energy  Storage  and  Attitude  Control 

The  basic  flywheel  energy  storage  and  attitude  control  unit  consists  of 
two  wheels  with  axes  carefully  aligned,  spun  in  opposite  directions  to  give 
the  desired  angular  momentum  and  energy  combination.  Note  that  only  the 
control  law  is  different  from  energy  storage  alone.  Associated  with  each  fly- 
wheel are  generating  and  motoring  electronics  and  suspension  electronics. 

System  control  electronics  are  common  to  both  wheels  and  create  wheel  speed 
change  commands  for  the  wheel  motor-generator  electronics. 

At  any  given  instant  of  time,  the  system  controller  desires  each  of  the 
two  wheels  to  be  running  at  a given  speed  to  contain  the  correct  total  energy 
and  angular  momentum.  Figure  1-1  gives  the  relationship  between  the  wheel 
speeds  and  o ^ anc*  the  angular  momentum  and  energy  of  the  system.  All  quan- 
tities are  normalized  to  the  maximum  speed  Note  that,  since  the  motor- 

generator  is  designed  to  operate  over  a limited  speed  range  (nominally  .A  u^<<jj<w^) 
only  the  upper  right  portion  of  the  plot  is  useful.  At  constant  stored  energy, 


TABLE  1-6 


REQUIREMENTS  FOR  LES  MOMENTUM  WHEEL  MAGNETIC  BEARINGS, 
ELECTRONICALLY  CIMBALLED  BASED  ON  TELDIX  MDR  WHEEL 
NOMINAL  MINIMUM  REQUIREMENT 


English 

MRS 

Angular  momentum 

> 75  ft. -lb. -sec. 

100  N-M-S 

Momentum  range 

> ± 10% 

± 10% 

Electronic  gimballing 

> 1.2° 

20  mrad 

Momentum  component 

> 1.4  ft. -lb. -sec 

2 N-M-S 

Slewing  rate 

> 5°/sec. 

0.1  rad/sec. 

Ginballing  torque 

> 3.5  ft. -lb. 

5 N-M 

Motor  torque 

> .035  ft. -lb. 

.05  N-M 

Nominal  speed 

> 16,000  RPM 

16,000  RPM 

Weight,  no  electronics 

< 22  lbs. 

10  kg 

OD 

< 12  inches 

. 3 m 

Height 

< 8 inches 

.2  m 

Power,  steady  state,  OG 

< 10  w 

10  W 

Caging  for  launch 

Yes 

Yes 

Suspension 

Electromagnetic,  5-axis 

active 

touchdown  ball  bearings 

Motor 

DC  brushless 

Redundancy 

All  motor  and  suspension  w 

indings . 

No  single-point  failures  except 
coil  shorts  or  mechanical  failure. 
Automatic  fault  detection  and  recovery. 
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the  system  operates  along  the  circular  arcs  centered  on  the  origin  to  vary 

angular  momentum.  Such  variations  wou Ld  be  quite  small  for  a system  with 

large  bias.  (Typically  H/Ito^  = .1,  Aw/uy.  < .01.)  At  constant  momentum,  the 

system  energy  is  charged  and  discharged  along  the  straight  lines  with  positive 

2 

slope  shown.  Note  that  for  these  values  with  H / T = 0,  AE/Iti ^ = .84  or  84% 
of  the  maximum  stored  energy  is  accessible.  For  H/Ioo^  = .1,  this  reduces  to  70%. 

It  is  comforting  to  know  that  such  a system  has  no  unusual  stability  prob- 
lems, and  indeed  the  control  system  shown  in  Fig.  1-4  represents  a first-cut  at 
a stable  control  system  diagram  allowing  energy  and  momentum  management  with 
relatively  simple  topology.  Such  a system  could  have  perhaps  10-sec.  transi- 
ent response  times  to  disturbances,  with  a second-order  response  and  adjustable 
damping  ratio.  In  reality  such  a system  must  be  carefully  modeled  to  include 
the  attitude  control  system  in  detail. 

Two  other  areas  of  concern  are  speed  control  and  wheel  alignment. 


Speed  Control 


Our  control  system  requires  that  the  speed  of  the  wheel  be  measured 
within  t = .5  sec.  Can  it  be  done  accurately  enough  in  that  time?  Assume 

(Jj 

the  motor-generator  has  N poles  which  are  sensed  as  they  pass.  Assume  a 
1 MHz  = f clock  counts  the  interval  between  poles,  so  that  in  t seconds  we 

U) 

receive  f t counts  and  the  wheel  rotates  N / 2tt  t poles,  00=271  (RPM)/60. 

Ob  00  0) 

We  can  then  measure  uj  within  one  part  in  fT  , which  for  our  numbers  Is 

5 * 8 

5 y 10  . This  is  well  within  the  stability  of  a good  clock  (parts  in  10  ) 

and  sets  the  angular  momentum  within  the  same  accuracy  for  constant  wheel 

inertia  I . An  error  of  H/fT  in  angular  momentum  causes  a satellite  rate 

“ 0)  6 3 2 3 

of  tij  = H/i  fT  which  for  our  numbers  is  I =10  Kg-m  , H = 10  N-M-S, 

s s Oi  s 

(iJfi  = 2 y 10  f’  which  is  well  below  the  operational  requirements.  Without 

feedback,  however,  knowledge  of  I within  this  accuracy  is  unlikely.  There- 
to 

fore,  if  the  ACS  loop  is  not  closed,  satellite  rates  of  oo  = HAI  /I  I may 

s w s 0) 

result . 
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Fig.  1-4.  Wheel  speed  control  loop  diagram 


-3  -3 

For  AI  /I  =10  , oj  = 10  or  .05° /sec.  may  result. 

0)  U)  a 

For  a turn  mode  at  .5°/sec.,  a 10%  error  will  result.  Stationary  constant 
speed  operation  for  several  minutes  would  result  in  several  degrees  drift  unless 
AI  were  measured  and  compensated  for. 

Al ignment 

If  the  two  wheels  of  a pair  are  misaligned  by  6 radians,  the  discharg- 
ing of  the  wheels  during  eclipse  could  transfer  .66  H ^eej_  angular  momentum  to 
the  satellite  over  1.2  hours  would  result  in  an  uncorrected  rate  of 

.66  H _ 

wheel 

Ws  (1.2)  (3600)  (H  bias) 

-A 

For  6 = .01°  = 2 * 10  rad,  which  appears  practical  for  a magnetic  bearing 

system,  we  get  with  H,  = 0.1  H . , 0)  = . 06°/1.2  hrs . = 1.4  x 10  ^ deg./ 

^ bias  wheel  s 

sec.  = 2. A x 10  rad. /sec.,  quite  low.  This  momentum  is  perpendicular  to 

the  base  momentum  of  the  satellite,  so  the  amplitude  of  the  resulting  motion 

is  limited  to  1.26  H . ,/M.  , . This  gives  .12°  maximum  motion,  which  is  quite 

wheel  bias  1 

acceptable  if  the  above  alignment  6 can  be  achieved.  Parameters  of  the  combined 
function  flywheel  system  are  shown  in  Table  1-7. 

The  results  of  the  system  trade-off  study  are  shown  in  Table  1-8.  The 
options  for  power  systems  are  batteries  I (NiCd,  no  discharger)  and  111  (Nill^, 
no  discharger)  from  the  power  study,  and  the  two  LES  flywheels  with  built-in 
momentum  storage.  Momentum  wheels  for  the  battery  systems  include  four  10  ft.- 
lb.-sec.  wheels,  one  100  ft. -lb. -sec.  wheel  without  redundancy,  one  10  ft.-lb.- 
sec.  wheel  plus  extra  gas,  and  a. 100  ft. -lb. -sec.  wheel  based  on  LES  technology. 

The  results  show,  as  expected,  that  all  NiCd  systems  are  very  heavy.  The 
difference  in  weight  between  the  lightest  N ill ^ system  and  the  lightest  LES-50 
flywheel  system  is  AA%,  which  shows  the  magnitude  of  payoff  in  technology 
investment  here.  The  LES-30  flywheel  is  29’/.  lighter  than  a NiH^  system  with 
conventional  wheel. 
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TABLE  1-7 


WEIGHT  AND  POWER  BUDGET  FOR  COMBINED  POWER  AND  ACS  FLYWHEELS 


Assume  2 Wheels  Total,  with  H'  - 

0.1: 

34  W-H/lb. 

54  W-H/lb. 

Characteristics  of  1 Wheel: 

O.D.  of  rotor.  In. 

20- in. 

20- in. 

RPM  Op.  range 

16K-4 OK 

20K-50K 

Rotor  Mass,  lb. 

20 

20 

Max.  stored  E,  watt-hrs.  (MJ) 

1360  (4.90) 

2160  (7.78) 

M/G  weight,  lb. 

10 

20 

Inertia,  Kg-m 

Available  H,  net  ang.  momentum. 

867  (121) 

1084  (152) 

ft  -lbs. -sec.  (N-M-sec.) 

89 

111 

Mag.  bearing  and  touchdown  wt. , lb. 

20 

20 

Vac.  housing,  lb. 

10 

10 

Shaft,  lb. 

6 

6 

Misc.,  lb. 

4 

4 

Rotating  unit  wt.,  lb. 

90 

100 

Available  E,  W-H 

952 

1512 

Available  E/M  W-H/lb. 

10.6 

15.1 

1.2  hr.  discharge  rate,  watts 

792 

1260 

Eclipse  Budget  per  2 Wheel  System 

Generator  efficiency  (Z) 

96 

96 

Bridge  efficiency  (Z) 

92.5 

92.5 

Drag  power  from  bus  (W) 

40 

40 

Suspension  and  control  power 

10 

10 

Output  power  to  bus  (W) 

653 

1069 

Generating  efficiency  (Z) 

84 

85 

Output  energy  to  load  (W-h) 

784 

1282 

Bridge  weight  .035  Ib./W 

21.4 

34.4 

Control  electric  weight  .005  lb./W 

3.0 

5.0 

Suspension  electric  weight  (lb.) 

5.0 

5.0 

Total  electric  weight  (lb.) 

29.4 

44.4 

Total  system  weight,  2 wheels 

119.4 

144.4 

Total  System  E/M  (W-H/lb.) 

6.6 

8.9 
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TABLE  1-7  (Cont'd) 


Sun  Load  Budget 

per 

2 Wheels-Charging , c/10  hr. 

Rate: 

34  W-H/lb . 

54  W-H/lb 

Maximum  drag  loss.es  (W) 

40 

40 

Suspension  and  contiol  (W) 

10 

10 

Wheel  Input  power  (W) 

00 

152 

Motoring  efficiency,  motor 

(7) 

945 

94.5 

Motoring  efficiency,  bridge 

(.7) 

89.5 

89.5 

Total  input  power  (W) 

173 

239 

Charging  efficiency  (%) 

55 

04 
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TABLE  1-8 


Note,  however,  that  there  are  two  ways  to  minimize  weight  once  the  NiCd 
battery  choice  is  eliminated.  One  can  go  to  a flywheel  storage  system  of 
advanced  design  for  the  lightest  option,  or  use  the  flywheel  technology  to 
build  a light  momentum  wheel,  shown  as  the  LL  100  FPS  wheel.  This  second 
option  is  far  easier  to  achieve,  and  gives  half  the  payoff  of  the  lighest 
choice.  The  key  parameter  becomes  high  wheel  speed,  resulting  in  low  weight. 
For  long  life,  this  requires  magnetic  bearings.  Such  systems  are  clearly 
available  in  the  time  frame  of  1980  (see  Teldix  MDR  momentum  wheel)  even  with 
steel  rotors.  The  addition  of  composite  rotors  then  makes  energy  storage 
feasible,  but  this  can  come  later. 

Integration  of  a flywheel  system  as  opposed  to  a battery  and  momentum 
wheel  combination  presents  advantages  for  each.  The  flywheel  approach  repre- 
sents a single  subsystem  to  be  built  or  purchased,  tested,  and  integrated 
instead  of  two  separate  ones,  representing  savings  in  cost  and  manpower  and 
perhaps  resulting  in  more  intense  refinement  for  the  same  resources.  However, 
power  and  attitude  control  now  interact  intimately.  Probably  two  sets  of 
flywheels  would  be  desirable,  one  for  ACS  system  tests  and  one  for  power  system 
simulation.  The  reliability  of  the  final  system  is  affected  both  ways  also  in 
that  there  is  only  one  subsystem  to  fail  rather  than  two.  Since  failure  of 
either  would  probably  end  the  mission,  this  is  an  improvement.  However,  making 
the  flywheel  reliable  depends  upon  the  solution  of  the  key  design  problems 
discussed  elsewhere. 

In  summary,  the  flywheel  energy  and  attitude  control  system  is  lighter 
than  any  other  currently  available  alternative.  It  matches  the  expected 
longevity  of  NiH^  batteries.  It  represents  a large  technology  investment  for 
relatively  small  weight  improvements  over  MiH^  battery  systems.  However,  much 
of  this  technology  should  be  developed  anyway  to  produce  lightweight,  magnetic- 
ally suspended  momentum  wheels.  Therefore,  the  technological  progress  of  these 
significantly  different  rotating  machines  can  be  phased  into  a lower-risk 
program  with  tangible  short-term  rewards. 

It  makes  sense  now  to  develop  a lightweight,  perhaps  electronically  gim- 
balled,  internally  redundant  high-speed  momentum  wheel.  With  that  in  hand, 
energy  storage  can  be  studied  later. 
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II.  TECHNOLOGY  SURVEY 


The  following  section  details  the  state-of-the-art  of  each  flywheel 
component . 

A.  Rotor 

For  any  given  system  parameters,  flywheel  system  energy  density  is 
roughly  proportional  to  that  of  the  rotor.  The  energy  per  unit  weight  which 
can  be  stored  in  a flywheel  rotor  is  directly  proportional  to  the  allowable 
stress  to  density  ratio  obtainable  in  the  material.  For  an  isotropic  constant 
strength  material,  such  as  some  metals,  the  mathematical  solution  to  a constant- 
stress  optimized  wheel  derived  by  Stodola  in  1924  still  represents  the  state- 
of-the-art.  This  gives  a wheel  with  thick  hub  with  an  exponentially  thinning 
rim,  which  is  truncated  in  various  ways  for  a practical  design.  However,  real 
materials  are  not  uniform  in  allowable  stress  because  their  properties  vary 
with  thickness  due  to  the  hardening  and  forming  processes  used.  Therefore,  a 
better  optimum  could  be  realized  by  including  this  dependence  in  the  mathe- 
matical optimization  of  the  shape.  However,  there  is  little  reason  to  believe 
this  optimum  will  significantly  improve  on  the  present  designs  (8  W-H/lb. 
obtained  with  titanium  Stodola  wheels).  One  such  wheel  system  with  mechanical 
bearings  was  developed  by  Rockwell  and  delivered  to  NASA/Langley  RC  this  year. 

Composite  filament  wound  materials  are  anistropic;  their  properties 
(strength,  stiffness,  etc.)  vary  with  orientation.  Tensile  strength  is  very 
high  in  the  direction  of  fiber  length  but  the  normal  strength  and  shear  strength 
are  quite  small.  Therefore,  it  is  possible  to  create  thin  hoop-shaped  struc- 
tures by  filament  or  tape  winding  processes  which  have  primarily  circumferential 
strength.  Thin  hoops  have  been  built  and  tested  which  demonstrated  80  to  90 
watt-hours  per  pound  using  Kevlar.  However,  attempts  to  make  the  hoop  radially 
thicker  produce  radial  stresses  in  the  weak  direction  and  failure  occurs  at 
lower  energy  densities.  This  is  often  aggravated  by  locked-in  thermal  radial 
stresses  due  to  processing  problems. 
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The  result  is  that  disc  or  hollow  cylinder  designs  tend  to  break  radially 
if  a significant  fraction  of  the  circumferential  design  strength  is  used.  j 

Simple  solid  circumferentially  wound  discs  of  Kevlar  composite  actually  have 
not  achieved  as  good  an  energy  density  as  isotropic  metal  wheels.  S-glass, 
which  has  somewhat  better  transverse  characteristics  due  to  a more  developed 
method  of  bonding  fiber  to  matrix,  does  somewhat  better  but  is  still  not 
impressive  for  solid  discs.  Coodyear  Aerospace  is  now  working  on  improved 
matrix  characteristics  for  Kevlar  wheels. 

One  approach  to  solving  this  problem  is  to  put  an  outer  radial  fiber  . 

wrap  of  Kevlar  on  the  rotor  after  the  circumferential  inner  hub  is  cured. 

This  adds  some  radial  strength  and  helps  carry  the  radial  loads  on  the  inner 

matrix  material.  Currently,  Rockwell  is  attempting  this  for  NASA/Langley , < 

so  far  without  success  beyond  5 watt-hours  per  pound.  At  present,  it  is  not 

clear  that  the  outer  wrap  does  not  create  more  problems  than  it  solves  by 

introducing  discontinuities  in  the  radial  properties. 

A second  approach  being  tried  is  to  build  up  a radially  thick  structure 
out  of  thin  hoops  which  are  discontinuous  and  thus  strain-relieved  at  their 
interfaces.  By  using  a number  of  concentric  thin  hoops  in  a bicycle-wheel 
shaped  assembly  with  flexible  spokes,  Carrett-Air esearch  Corporation  has 
achieved  an  impressive  17-watt  hour  per  pound  performance  on  a number  of 
Kevlar  rotors  storing  from  tens  of  watt-hours  to  several  kilowatt  hours. 

The  precise  method  of  bonding  concentric  rings  together  at  Garrett  is  propri- 
etary but  involving  preloading  fiber  tension  and  ring  interference  fits  held 
together  by  friction.  However,  these  structures  were  test  devices  and  need 
to  be  evaluated  for  flight  criteria.  A variation  of  this  approach  using 
concentric  rings  of  different  properties  tied  together  with  elastomers  has 
been  plagued  with  dynamics  problems.  Dynamics  problems  broke  the  wheel-hub 
bond  in  the  Garrett  wheel. 

A third  approach  is  to  weave  the  fibers  in  the  rotor  in  different  direc-  i 

tions  so  as  to  give  the  structure  strength  in  the  radial  and  axial  directions 

I 

as  well  as  the  circumferential.  This  requires  a 3-dimensional  "polar  weave" 
capability  which  General  Electric '6  RESD  facility  and  others  have  developed 
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for  nose  cones.  CE  has  produced  several  Kevlar  rotors  of  this  design  and 
expects  to  achieve  40  watt-hours  per  pound  with  them,  but  have  not  yet  tested 
them  beyond  the  b watt-hour  per  pound  level.  Testing  to  the  higher  level  was 
scheduled  for  December,  1976.  This  may  also  be  a viable  solution,  but  is  very 
expensive . 

A fourth  variation,  used  by  AVCO  and  developed  to  about  the  same  point 
over  a longer  period  of  time,  builds  up  a radially  thick  flywheel  by  using  a 
two-dimensional  wea\e  to  get  circumferential  and  radial  strength  and  then  lam- 
inate these  together  axially.  These  separated  from  the  hub  at  8 watt-hours/lb. 

A3  1 four  of  these  approaches  are  currently  under  test  and  should  show 
20  to  40  watt-hour  per  pound  performance  if  design  expectations  are  achieved. 
Many  of  the  current  tests  are  motivated  by  U.  S.  Army  Ft.  Belvoir  contracts  for 
mobile  energy  storage  devices  under  the  MERADCOM  program. 

The  most  promising  current  experiments  with  anistropic  rotor  design  are 
being  done  at  Johns  Hopkins  Applied  Physics  Laboratory.  Using  bare  filament 
Kevlar  and  other  materials,  rotors  are  wound  without  glue  on  a metal  hub. 

Spokes  of  the  same  material  are  wrapped  through  the  hub  and  around  the  peri- 
phery, then  glued  to  give  some  radial  and  axial  strength  to  the  wheel.  The 
resulting  rotor  is  attached  to  a thin  shaft  with  a metal  and  rubber  shock 
mount  to  allow  self-centering.  Results  have  reached  30  to  40  watt-hours  per 
pound  in  this  solid  disc  configuration,  in  spite  of  problems  in  bonding  spokes 
to  circumferent ials.  Due  to  the  lack  of  glue  on  the  circumferentials , the 
potential  hoop  strength  is  nearly  double  other  designs. 

In  summary,  the  composite  rotor  represents  a technical  risk  but  appears 
well  worthwhile  as  a weight-saving  advantage  and  an  advance  in  the  state-of- 
the-art  for  energy  storage. 

B . Magnetic  Bearings 

The  magnetic  bearing  represents  a quantum  improvement  in  the  long-term 
reliability  and  lifetime  expectancy  of  high-speed  spinning  assemblies.  A 
combination  of  permanent  magnetics  and  electromagnets  is  used  to  levitate  a 
rotor  between  the  poles  of  the  magn<  tic  bearing.  The  resulting  structure  has 
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no  contacting  surfaces  to  require  lubrication.  The  bearing  requires  electrical 
power  only  to  counteract  disturbances  and  can  be  operated  with  no  steady-state 
dissipation  at  its  bias  point.  The  bearing  is  limited  in  force  density  by  the 
magnetic  saturation  characteristics  of  iron,  so  that  size,  weight  and  ground 
power  are  dominated  by  the  need  to  support  the  rotor  during  testing  in  a 1G 
env  ironment . 

In  addition,  a magnetic  bearing  is  a "soft"  support  in  that  its  restoring 
spring  constant  is  relatively  weak  compared  to  mechanical  bearings.  Therefore, 
such  a system  is  characterized  by  clearance  dimensions  measured  in  mils  or 
tens  of  mils,  rather  than  the  tens  of  micro-inches  typical  of  conventional 
ball  bearings.  Such  soft  bearings  effectively  isolate  rotor  from  mounting 
when  imbalance  vibrations  occur.  Also,  they  allow  a relatively  crudely  bal- 
anced rotor  to  spin  about  its  principal  axis  with  minimum  effect  on  the  system 
performance.  However,  the  low  stiffness  places  limits  on  the  allowable  cross- 
axis spin  rate  which  causes  lateral  torque  on  the  spinning  flywheel. 

At  present,  several  gyro  and  momentum  wheel  vendors  are  investing  in 
magnetic  bearing  technology  to  apply  these  bearings  to  inertial  navigational 
systems  and  large  momentum  wheels.  Cambion  Corporation,  which  invented  the 
magnetic  bearing  and  has  accumulated  a long  record  of  innovations  in  this  field, 
is  retrofitting  several  wheels  for  Bendix  in  hopes  of  interesting  some  customers 
of  Bendix  systems  in  an  experiment  promising  longer  lifetimes  and  eventually 
higher  speeds.  Sperry-Rand  Corporation,  Flight  Systems  Division,  has  developed 
a considerable  in-house  magnetic  bearing  capability  and  is  rapidly  developing 
these  bearings  for  its  momentum  wheels  to  extend  their  lifetime. 

Most  work  at  Sperry  has  been  with  radially  passive,  axially  active  perman- 
ent magnet  systems  for  momentum  wheels.  However,  the  wide  servo  bandwidths  and 
associated  stability  problems  of  these  systems  have  caused  recent  work  to  be 
directed  at  5-axis  active  electromagnet  systems. 

Aerospatiale  (France)  has  developed  a permanent  magnet  bearing  for  a 
momentum  wheel  with  a composite  rotor.  This  radially  passive,  axially  active 
system  has  realtively  high  input  power  requirements,  due  in  part  to  the  imbal- 
ances from  a composite  rotor  and  in  part  from  the  many  magnetic  components  in 
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this  design.  It  includes  a passive  magnetic  damper  for  oscillations  about  the 
transverse  directions.  This  wheel  design  is  an  ambitious  use  of  several  new 
technologies,  with  a rather  complex  implementation.  It  is  the  most  developed 
unit  of  its  kind,  and  has  been  qualified  for  COMSAT. 

Teldix  (Germany)  has  a metal  wheel  with  5-axis  active  magnetic  bearings. 
This  allows  limited  magnetic  gimballing  of  the  rotor  by  changing  the  set  points 
of  the  transverse  position  loops.  The  limit  is  the  maximum  gimbal  angle,  done 
to  small  clearance  in  the  magnetic  gaps.  The  KDR  series  wheel  at  100  N-M-5  has 
a ± 1/2°  gimbal  capability,  and  has  many  coils  for  sensors,  bearings,  and  motor 
which  are  single-point  failures.  However,  engineering  radial  wheels  with  1.2° 
gimballing  and  fully  redundant  coils  have  been  built  (not  qualified)  and  demon- 
strated automatic  fault  correction.  The  MDR  wheel,  flight  qualified  for  COMSAT, 
appears  to  be  a very  well-engineered  unit,  clearly  implemented  and  easily 
inspected  or  repaired.  The  bandwidth  of  its  servo  loops  is  set  by  the  wheel 
nutation  frequency,  and  is  not  too  much  lower  than  the  comparable  permanent 
magnet  system  due  to  high  (16K  rpm)  wheel  speed.  Other  manufacturers  of  high- 
speed rotating  equipment  are  moving  more  slowly  in  applying  these  bearings. 

A fundamental  statement  derived  from  Maxwell's  equations,  known  as 
Earnshaw's  Law,  states  that  no  static  assembly  of  permanent  magnets  and  ferro- 
magnetic materials  will  suspend  a device  in  stable  equilibrium.  This  means 
that  at  least  one  degree  of  freedom  of  the  rotor  in  a magnetic  bearing  system 
must  be  stabilized  by  active  feedback  electronics.  Passively  stable  behavior 
can  be  obtained  axially  (along  the  spin  axis)  or  radially  (perpendicular  to 
the  spin  axis)  but  not  both.  For  energy  storage,  it  appears  that  a passive- 
radial,  active-axial  system  will  allow  fully  redundant  electronics  and 
electromagnets  to  be  used  with  a simple  arrangement  described  in  Section  4. 

The  radially  passive  containment  can  be  achieved  with  permanent  magnets 
while  electromagnets  at  each  end  of  the  shaft  are  used  for  axial  stabiliza- 
tion. Proper  design  allows  either  electromagnet  to  perform  the  function  if 
one  should  fail.  Key  design  parameters  are  low  weight  and  realization  of 
the  potential  for  low  power.  This  design  represents  state-of-the-art  engin- 
eering using  Samarium-Cobalt  permanent  magnet  materials  in  ways  which  do 
not  exceed  their  limited  stress  capabilities. 
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C . Touchdown  Bearings 

The  magnetic  bearings  would  not  be  designed  to  withstand  launch  vibration 
and  shock,  power  failures,  or  accidental  test  handling  mishaps.  These  loads 
would  be  absorbed  by  a "touchdown"  or  backup  ball  bearing  which  makes  no  con- 
tact with  the  rotor  under  operating  conditions,  but  prevents  contact  of  the 
magnetic  components  and  complete  failure  in  case  of  mishap.  This  bearing  will 
rerjuire  careful  design  to  avoid  brinelling  or  breaking  during  a high-speed  spin 
down.  It  must  withstand  high  accelerations  without  significant  deterioration. 

For  large  energy  wheels,  this  bearing  represents  a development  item  which 
Sperry,  Bendix,  or  Draper  Laboratories  can  design  with  existing  technology  and 
is  an  engineering  development  item  of  relatively  low  risk.  For  smaller  wheels, 

the  present  MIT/LL  design  is  adequate.  j 

i 

D . Mo tor -Genera tor 

The  design  of  small,  light,  efficient  motor  generators  is  an  art.  This 
art  is  understood  by  a number  of  manufacturers  including  GE  and  Garrett.  The 
key  to  this  design  is  the  use  of  Samarium-Cobalt  magnets  on  the  rotor  without 
exceeding  the  stress  limits  on  the  material,  and  the  minimization  of  eddy  cur- 
rent losses  in  the  stator  windings  without  sacrificing  efficiency.  Seventy 
percent  cycle  efficiency  for  a 10-hour  charge,  12.8  hour  coast,  1.2  hour  dis- 
charge is  desired.  Rockwell's  flywheel  uses  a GE  design  and  achieves  adequate 
motor-generator  performance,  so  there  is  no  reason  to  expect  a technical  risk 
in  this  item.  The  current  effort  with  Sperry  is  reported  elsewhere,  and  is  a 
better  solution  to  the  problem. 

E . Shaft  and  Housing  Design 

The  shaft  and  housing  for  the  flywheel  system  represent  a significant 
part  of  the  system  weight.  They  are  to  provide  rigid  mounting  support  for 
wheel  alignment,  estimated  to  be  tens  of  arc  seconds  to  avoid  attitude  control 
interaction  with  power  demands,  and  also  provide  a vacuum  enclosure  for  ground 
testing.  There  is  no  doubt  that  these  can  be  built.  Use  of  composite  materials 
for  lighter  weight  and  perhaps  flywheel  failure  containment  will  be  considered 
but  may  be  a technical  risk. 
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III.  ROTOR  ANALYSIS 


A.  Orthotropic  Rotor , Analytical  Analysis 


Equilibrium  equation: 


(hra  ) - haQ  + pu)2hr2  = 0 
dr  r o 


Compatibility  equation: 
deQ 

r -j — + eQ  - e =0 
dr  9 r 


Fig.  3-1.  Orthotropic  rotor  geometry. 

Stress-strain  (plane  stress) : 

Cr  IT  ' VrV’  £0  = e7  (°6  ' Vr> 

r 0 

Introduce  the  stress  function,  tjr. 


_ , . dtp  2.  2 

ha  - ip,  haQ  = + pa)  hr 

r r 0 dr 


and  eliminate  e , eQ  from  the  compatibility  equation,  by  means  of  the  stress- 
r o 

strain  relations: 


d r 1 ,1  2 2 tj;  . , 
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1 , 1 d^  2 2 ip  . 
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arc  constants: 


If  E0>  £r’  V \ 


rV'  ♦ rip'  fl  - vf)  + Y~  vr  " TT  r] 


r f. 

* <P  lvB  E“  r " 1 = -(7,  + T-  V pU)2r^h 

r Jr 


In  the  general  case  (neglect  variation  of  v„) 


1-V'  + [1  ■ r^T‘  + ^ Vr  ■ fTrl 


F " F. 

, f , h'  0,  0 

^(ve(rr+rr)T 


^ fl  2 3 

= .(J-r^^vr)  Pu>  r h 

^0  r 


y,  y2  . 

NOTH:  If  h “ r , E0  « r , eQ/er  const,  the  equation  remains 


CASF.  1:  h = const.  (=1),  eQ,  cr,  vfi,  vr 


const . 


r2<p"  + r ip'  [1  - v0  + jr1  V - pT  = “(S  + V p°rr 
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, 2o  Eo 

Let  Vr  I"  = V E"  = n 
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2 3 


r2\p"  * r ip'  [1  - VQ  + V ] - i^n  = - (3  + vff)  pto  r 


equidimensional  . 


33 


ip"  + V i\  - v0)  - = -(3  + \*r)  pu  e' 


y...» 


General  notation: 
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B . Bare  Filament  Rotor  Analysis 
1.  Basic  Assumptions 

The  configuration  of  a bare  Kevlar  fiber  flywheel  under  steady  rota- 
tion is  determined  in  this  section.  It  is  supposed  that  the  evenly-spaced 
radial  spokes  are  inextensible , although  the  following  analysis  may  be  readily 
extended  to  Include  spoke  compliance.  Furthermore,  the  individual  fibers  compos- 
ing the  circumferential  windings  are  assumed  to  be  mutually  noninteracting 
(i.e.,  no  contact  or  friction  forces  are  generated  between  fibers)  so  that  they 
may  be  considered  separately.  Finally,  each  fiber  is  approximately  a homogene- 
ous elastic  string  with  zero  torsional  and  bending  rigidity. 

In  the  next  section,  the  equations  for  the  steady  rotary  motion  of  a string 
are  assembled  and  the  boundary  conditions  for  the  bare  fiber  flywheel  are  derived. 
Next,  the  approximation  which  neglects  the  compliance  of  the  fibers  is  justified, 
and  explicit  results  are  obtained  for  the  configuration  and  tension  distribution. 
After  presenting  certain  useful  approximations,  we  conclude  with  a calculation 
of  the  loads  Imposed  on  the  spokes  by  circumferential  windings. 
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2.  Equations  of  Steady  Rotary  Motion 


Following  Ref.  (1),  the  general  equations  of  motion  of  a homogeneous 

\ 

elastic  string  may  be  listed  as  follows: 


JL  / r i2L 

as  \ 9s 


/a  . a \ /ax1  ax1\  . _i 

m (at  + u [W  + u aT  + F 


(3-1. a) 


Tangent  Identity 


ax1  3X1  = 1 

at  3s 


(3-1. b) 


Mass  Conservation  Equation 


am  , 3(mu)  _ n 
at  a8 


(3-1. c) 


Stress-Strain  Relation 


m - mQ  exp(-kT) 


(3-1. d) 


where  we  define: 


s ■ arclength  of  stretched  string  (measured  from  an 
attacliment  point) 

£ 

sq  » arclength  of  unstretched  string 
i A 

x « inertial  position  vector  of  a string  point 

u - tangential  component  of  velocity  of  string  point  at 

fixed  s relative  to  a point  at  fixed  s. 
o 

m * lineal  mass  density  of  stretched  string 

k 

- ra(T-O)  * lineal  denstiy  of  unstretched  string 


T - tension 


k ■ axial  compliance  of  string  materia] 
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and  where  F*  is  the  externally  applied  force  per  unit  length.  Since  the  Inter- 
action of  each  fiber  with  other  fibers  is  neglected,  F*  - 0 in  this  analysis. 
Also  note  that  the  stress-strain  relation  implies  the  use  of  the  logarithmic 
or  natural  strain  and  is  valid  for  finite  strain. 

Consider  the  steady  rotary  motion  of  the  flywheel  with  angular  speed  0 
about  the  flywheel  axis.  In  this  case,  each  fiber  will  lie  in  a plane  normal 
to  the  rotation  axis.  Because  the  motion  is  steady,  u vanishes  and,  according 
to  the  mass  conservation  relation,  m is  independent  of  time.  The  acceleration 
of  a string  point  thus  becomes: 


3 


2 i 
x 


-a 


2 i 
x 


Adopting  polar  coordinates  (r,  41)  in  a plane  rotating  with  the  flywheel,  the 
equations  of  motion  become  (see  Reference  (2)  for  comparison). 

Radial  Momentum: 

(Tr')’  - Tr<f«'2  + mnZr  = 0 (3-2. a) 


Circumferential  Momentum: 

[Tr24>']  - 0 

Tangential  Momentum: 

T'  + mii2rr'  » 0 

Tangent  Identity: 

r’2  + r <J>  ’ 2 = 1 

Stress-Strain  Relation: 

m = m exp  (-kT) 
o 

Also,  from  mass  conservation,  we  have: 
m 

0 8 O 

— - — 

r)s  m 

o 


(3-2. b) 

(3-2. c) 

(3-2. d) 

(3-2. e) 

(3-2. f) 
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Finally,  it  is  most  convenient  to  adopt  <\>,  rather  than  s,  as  the  independent 
variable.  Using  the  tangent  identity: 


r2  4 


the  momentum  equatiors  become: 


<L_  Tr2 

d«J> 


r2  + 


(3-3. a) 


dT  n2  dr 

^ + mT2  r — - 0 


(3-3. b) 


where  the  radial  momentum  equation  may  be  omitted  since  it  is  merely  a combina- 
tion of  the  circumferential  and  tangential  equations.  Note  also  that  the  element 
of  unstretched  length,  ds  , ia  obtained  by  combining  (2.f)  with  the  relation 

for  M • 


, m + 2 

ds  =*  — wr  + 
o m i 
o 


(3-3. c) 


Relations  (3),  along  with  the  stress-strain  relation,  (2.e),  constitute  a 
complete  set  of  equations  for  each  fiber. 

To  derive  the  appropriate  boundary  conditions,  consider  a single  fiber 
attached  at  radial  distance  r^  along  the  neighboring  spokes.  Referring  to 
Fig.  3-2,  we  define: 

$ “ angle  between  spokes  (a  submultiple  of  2ir) 

L “ unstretched  length  between  spokes 
Because  of  symmetry  of  the  fiber  configuration  about  <}>  = <!>/ 2 : 
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I 


4)  = 0 : 

i p » <t/ 2 : 


r « r 


dr 


d <t> 


(3-4. a) 
(3-4. b) 


In  addition,  the  total  unstretched  length  must  be  constrained  to  be  L. 
Prom  (3.c): 


M2 

f 


d«t 


^Vr2  + 

m 


) -i 


(3-4. c) 


Equations  (4)  now  constitute  the  necessary  boundary  conditions. 

Before  integrating  Eqs.  (3)  and  (2.e),  it  is  convenient  to  define  the 
dimensionless  variables  T,  r,  v,  l,  and  K by: 


n2  2 
m Si  r 
o o 


r = r/r 

o 

~2  (3-5) 

v = r 

St  - L/r  4> 
o 

. n2  2 
k = km  r 
o o 

Clearly,  T and  K are  the  dimensionless  tension  and  compliance  while  9.  is  the 
ratio  of  the  total  length  to  the  arclength  of  a circular  arc  between  the  spokes. 
It  is  expected  that  as  a result  of  the  fibrication  process,  there  will  be  excess 
fiber  so  that  St-1  > 0. 

In  terms  of  the  dimensionless  variables,  the  governing  equations  and 
boundary  conditions  become: 


d_ 

d<J> 


TV 


VA1_  dv 

v 4v  d<)> 


i-  0 


(3-6. a) 
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'.v  ' 


d / 1 IfT 
d<t>  \ K G 

♦M 

- 0 

(3-6. b) 

= o 

: v = 

1 

(3-6. c) 

. * 

<P  - 2 

dv 

: d<}>  ~ 

0 

(3-6. d) 

$/2 

1 

f 2 

| d<t> 

*/ 

-KT  \ 

e \ 

, 1_  /dv  \ 
v 4v  \d<}>/ 

M/2 

(3-6. e) 

o 

Solution  of  the  Inextensible  Case 

Note  that 

the  term 

<T  is  the  axial 

strain. 

Supposing  that  TK  < M, 

M > 0,  Kqs.  6(b)  and  6(c)  become,  upon  rearrangement: 


d_ 
d ip 


t(1  + 


V 


0 


V2  I 2~ 

if  (SO  ' I(1  + 02> 

o * 

where  and  a 2 are  both  positive  and: 

«!  < 5 («M  - 1 - H) 
a2  < M 

Since  the  strain  at  failure  for  KEVLAR  29  is  less  than  four  percent,  M < . 04 
and 

< .02 

«2  < . 04 
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Thus,  as  a first  approximation,  terms  of  the  order  of  the  compliance  may  be 
neglected  and  the  fiber  treated  as  inextensible . 

With  this  approximation,  Eqs.  6(a)  and  6(b)  may  be  integrated  to  give: 

T = T + ^ (1-v)  (a) 

O Z 

(3-7) 

— = v Cos  ip  (b) 

r ° 

v + 1- I**) 

4v  v dip  ) 

where: 

Tq  ” T((f  - 0) 

ij,  ” angle  between  the  tangent  and  the  <P  direction  at  <}>  = 0. 
Elimination  of  t from  Eq.  (7)  gives  the  equation  determining  v(<t): 


f— ) 


2 2. 
T Cos  ip 
o 


F(v) 


4 


F (v)  » | 4 v(Tq  + (1-v))2  - To2Q>s“i|;| 


(a) 


(b) 


(3-8) 


- (v  - Vj)  (v2  - v)  (v3  - v) 


V1  - V2  - v3 


The  roots,  v. , v2>  and  v^  of  F(v)  determine  the  extrema  of  the  string 
curve  and  may  be  evaluated  iteratively  from  the  expressions 


2t 

vl,2  “ 1 + At  -1 


T -1  +V  (T  -l)2  + 

O 10 


(At  -1) 
o 


l«*  ♦ i M 1 

T 

o 


T (T  -1)  T 2Sin2^ 

„ , , o o . o 

V - 4t  - 4 r 4 - - „ 

3 ° ’s'1  (»3  - l)2 


ka) 


(b) 


(3-9) 
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to  yield  the  numerical  results  shown  in  Fig.  3-3. 

It  is  seen  from  8(a)  that  is  defined  only  for  v £ [v,,v~]  and  v > v_. 

d<p  12  — 3 

Since  the  latter  possibility  corresponds  to  an  unbounded  string  curve,  it  must 
be  discarded.  Furthermore,  Fig.  3-3  shows  that  < 1 and  v^  > 1.  Thus,  in  view 
of  the  boundary  conditions  6(c)  and  6(d),  the  string  curve  must  be  confined  to 
v £ [l,v  ],  and  6(d)  merely  implies  that  v(<J)  * i>/2)  = v 

Integrating  8(a)  and  expressing  $ as  a function  of  v,  we  have 


<J)  = TqCos^ 


{ x^F(x) 


, or  (Ref.  3,  p.  227) : 


2toCos^ 


V1  \V3_V1 


II  sin 


rvi  1 _ v_i  /vvi 

V2'V1  ’ V1  ’\V3_V1 


I \ v2_vl  V1  \v3  V1 


(3-10) 


where  II (<t> , n,  k)  is  the  elliptic  integral  of  the  third  kind.  The  condition 
v(<(  = Q/2)  = v 2 gives: 


^ 2toCosi|; 


— n I ! . _2 
I l 2 ’ v.  ’ 

V1  VV3~V1  \ 1 


IlZi 

lv3-vl 


IT  sin 


v2_vl 


, 1 - 


(3-11) 


This  provides  one  relation  with  which  to  determine  tq  and  \l>.  The  second 
relation  is  Eq.  6(e)  which,  when  the  compliance  is  neglected,  becomes: 
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f / 


, T + ~ (1-V) 

dv  o 2 


^(v3-v)  (v2-v)  (v-Vj  )‘ 


Or,  (Ref.  3,  pp.  219-220): 


(tq  + \ (i-v^))  F(6,  q)  - |'(v2~v3)  n(<5>  92,  q) 


■ i M ^7 


(3-12) 


where : 


A , -l 

~ sin  4 


^V3"V1)(V2_1) 


(V2— V1 ) (v3_1) 


q 


i_ 

l 


and  where  F(d,  q)  is  the  elliptic  integral  of  the  first  kind. 

The  quantities  and  determined  by  (11)  and  (12),  are  of  importance 
since  is  the  maximum  tension  in  the  string  (see  F,q.  7(a))  while  the  radial 
load  applies  to  the  spokes,  T^siniJ;,  is  proportional  to  both  and  tp. 

4.  Approximate  Evaluation  of  and  tj; 

To  obtain  useful  approximate  results,  note  that  when  9.  is  ^ 1, 
is  also  ^ 1,  and  the  factor  v^-x  in  the  integrands  leading  to  (11)  and  (12) 
can  be  approximated  by  (v^-1).  Thus,  for  | J?,—  1 1 « 1,  the  conditions 
v(<p  = $/2)  * v 2 and  6(e)  become  approximately: 
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1 1 T«cos^  / 


J x ^(v2-x)  (x-vx) 


/■2  T + 7 (1-x)  

2 / dx  ■■  — — * 4*  J^l 

-/  ^(v2~x)  (x-v^ 


Or,  evaluating  the  integrals: 


T Cosijj 


/Vlv2(v3- 1) 


tt  . -1  / V1+V2_2v1V2 

tt  - sin  I 

2 \ v2-v2 


. / . i/.  VMW  . , -i  /“i^-2 
2 \ To  + T \2  2 ))  T + Bin  h^- 


^”3-1  + ^(v2-l)(l-v1) 


(3-13) 


These  relations  may  be  further  simplified  by  approximately  v2,  and  v3. 
With  | 4-1 | « 1,  both  | Tq— 1 | and  ^ are  also  small,  and  Eq.  (9)  give  as  a first 
approximation: 


V1  1 + A " £ 


V2  ^ 1 + A + £ 


where 


A - j (Vl) 


( 3-14) 


| ^(xo-l)2  + 3S.  2i]> 


v3  2^  4 


x 


Expressing  t and  ij)  as  well  as  v1  and  v0  in  terms  of  A and  e and  substituting 
(14)  into  Eq.  (13),  we  finally  obtain: 


T s,  1 + -2  A 

o — 2 

Sim l^yJJ  ^C2  - A2 
with  A and  e determined  by 


(a) 

(b) 


e 


crc 


1-c. 


2 

2 


(e) 


(3-15) 


A - -C^c 


(d) 


where 


cos  y 


4 


2 2 

3 (1+A)  - e 


(1+|A)  ^l-|(e2-A2) 


(e) 


cos 


V~3  + Ve 


- A 


2(1  + A) 


(f) 


Numerical  results  for  T and  \p  as  functions  of  the  excess  length  (0.-1)  are 

7T  2 TT 

shown  in  Fig.  3-4  for  4,  5,  and  6 spokes  (i.e.,  4>  - y,  yir,  y radians,  resp.). 
Although  the  accuracy  degrades  somewhat  above  (£-1)  ^ .2,  these  curves  still 
reflect  the  correct  qualitative  trends.  It  is  seen  that  the  attachment  angle, 
ip,  is  roughly  the  Bame  for  all  values  of  4>  considered.  As  should  be  expected, 

the  maximum  tension  i falls  below  1 for  small  (0-1)  and  ultimately  rises 

o 

above  this  value  for  sufficiently  large  (0-1).  Although  the  maximum  is 
seen  to  decrease  somewhat  with  increasing  number  of  spokeB,  t n,  1 (which 
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0*  I* 


0 1 0.2  0 3 0«  0 5 


(/-I) 

Fig.  3-4.  Maximum  tension  and  attachment  angle. 
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corresponds  to  the  "hoop"  stress  of  a circular  arc  of  string)  is  still  a 
reasonable  approximation  over  the  range  of  parameters  considered. 

5.  Computation  of  Spoke  Loads 

Before  calculating  the  loads  applied  by  the  fibers  to  the  spokes  as  a 
result  of  Eq.  (15),  we  must  first  determine  the  unstretched  length  of  all  fiber 
segments  between  the  spokes.  When  the  spokes  are  completed  and  the  end-plates 
used  during  fabrication  are  removed,  the  lateral  surfaces  of  the  spokes  must 
support  a load  due  to  the  circumferential  windings.  This  lateral  load  causes 
the  spokes  to  thLcken  along  the  flywheel  axis  with  a corresponding  decrease  in 
the  radial  extent  of  the  spokes.  It  is  reasonable  to  assume  that  the  radial 
position  of  attachment  of  each  fiber  segment  decreases  proportionally,  '"hat 
is,  if: 


r (n)  ~ radius  of  n^  winding  during  fabrication 
o 

r (n)  ^ final  radius  of  n^  winding  (at  57.  = 0) 
o 

than  our  assumption  is  that  — (r  (n)  - r (n))  is  the  same  for  all  windings. 

r (n)  ° ° 

o 

Thus,  the  unstretched  length,  L(n),  of  a fiber  segment  belonging  to  the  nt*1 
winding  is  <t>ro(n),  and  the  dimensionless  quantity  becomes: 

f.(n)  * £,  independent  of  n. 

With  this  assumption,  both  the  dimensionless  tension  T and  the  attachment 
angle  are  the  same  for  all  fiber  segments.  The  tension  in  a fiber  at  its 
point  of  attachment  to  a spoke  is  then: 

m 572 r 2t  (d,  <&) 
o o o 

so  that  the  total  radial  force  on  the  spoke  due  to  one  fiber  Is: 

2SiniJmo^2ro2To(«,>  4>) 
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so  that  the  total  radial  force  on  the  spoke  due  to  one  fiber  is: 


2Sin4rniQn^ro2To(J,,  4>) 


If  the  circumferential  fibers  are  held  in  a hexagonal  close-packed  pattern 

by  the  spokes,  the  number  of  fibers  per  unit  cross-sectional  area  of  the  fiber 
2 1 

pack  is  r — 2",  where  d,.  is  the  fiber  diameter. 

\J3  df 


Consequently,  the  total  radial  load,  F^,  applied  to  a spoke  per  unit  radial 
distance  is: 


——  -^r  m S72r  ^ S . tpT  (H,  <£) 
, 2 o o o 


JT  d 


(3-16) 


where: 

rQ  “ radial  coordinate  along  spoke 
hw  “ thickness  of  flywheel 
dj.  ” fiber  diameter 

The  total  radially  directed  force  resultant,  T , developed  over  spoke 
cross-sections  as  a result  of  the  above  load  distribution  may  be  determined  from: 


— + F + M f!2r  - 0 

: r so 


V fRH*  V 


(3-17) 


T (r  - R ) - 0 
s o w 


where: 


T — force  resultant  over  spoke  cross-section 

8 

“ mass  of  spoke  per  unit  length 
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R “ 


RH 


flywheel  radius 

radial  coordinate  of  hub-spoke  attachment. 


Substituting  relation  (16)  for  F , Eq.  (17)  yields: 


T (r  ) 
s o 


1 2 2 
4 M JT(R 
2 s w 


r : 2) 


+ — - — rn  J12(R  3 * - r 3)  SinipT  (Z,  4>) 

JyfT  d/  0 ” 


(3-18) 


Thus,  the  total  load  applied  to  the  hub  by  one  spoke  is  Tg(ro  = R^) • 

The  first  term  in  (18)  results  from  the  certrifugal  loading  only  and  the 

corresponding  tension  in  individual  spoke  fibers  is  no  more  than  half  the  "hoop" 

2 2 

tension,  m^ft  R^  , developed  in  the  outer  circumferential  fibers.  Thus,  it  is 
the  second  term  in  (18),  representing  circumferential  fiber  loading,  that 
significantly  degrades  the  maximum  energy  density  below  the  optimum.  Denoting 
by  T^  the  contribution  of  this  term  to  the  total  load  applied  to  the  hub,  we 
have: 


h - rj=  - V>  S‘"*To(*>  « (3-») 

3V3  df 

Assuming  the  following  typical  values: 

R =9  in. , R - 2 in. 

W H 

h * 4 in. , d.  « 15  mil 

w 2 f 

3 

fiber  density  » .052  lbm/in 

2 

Figure  3-5  shows  T /fl  as  a function  of  (£-1)  for  A,  5,  and  6 spokes.  It  is 
seen  that  while  some  slight  advantage  is  derived  by  increasing  the  number  of 

spokes,  the  only  effective  way  to  eliminate  T^  is  to  reduce  the  excess  length, 

(Z-l),  as  far  as  possible. 
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(/-I) 

Fig.  3-5.  Contribution  of  circumferential  fiber  loads  to  total 
load  applied  to  the  hub. 


6.  Conclusions 

The  above  analysis  shows  that  both  the  fiber  configuration  and  tension 

depend  primarily  on  the  "excess  length"  of  the  fiber  segments  (i.e.,  the  difference 

between  the  true  length  and  the  length  of  a circular  arc  joining  the  fiber-spoke 

attachment  points) . While  the  maximum  tension  decreases  somewhat  with  increasing 

2 2 

number  of  spokes,  it  is  still  of  the  order  of  the  "hoop"  tension,  mQJl  rQ  . Fur- 
thermore, the  angle  formed  between  a fiber  and  the  circumferential  direction  at  a 
spoke  is  nearly  independent  of  the  number  of  spokes  and  linearly  dependent  on  the 
excess  length. 

Since  the  excess  length  ultimately  produces  significant  loading  of  the  spokes, 
flywheel  optimization  depends  on  the  maximum  reduction  of  this  quantity. 
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(V.  MAGNEifC  BKAR  INC.  INVESTIGATION 


A design  layout  of  the  MIT/LL  magnetic  bearing  test  bed  system  is  shown 
in  Tig.  4-1.  It  consists  primarily  of  magnetic  bearings  (stators  and  rotors), 
touchdown  bearings,  motor-generator,  turbi ne- f ] ywheel , shaft,  position  sensor, 
and  test  bed  support. 

The  following  sections  are  the  supporting  design  calculations  for  the 
magnetic  bearing  brassboard  and  descriptions  of  other  related  work. 

A.  Meehan i cal  Design  Goncept 

The  prototype  energy  storage  wheel,  assembly  consists  of  rotor,  shaft,  two 
magnetic  bearing  assemblies,  two  touch-down  ball  bearing  assemblies,  a motor 
generator  assembly  and  a base  support  structure.  A conventional  rotating  shaft 
spindle  configuration  vs.  fixed  shaft  was  chosen  for  the  prototype  model  as 
shown  in  Fig.  4-1.  The  rotor  was  designed  to  simulate  the  mass  of  a flight  con- 
figured wheel,  but  not  the  inertia.  The  rotor  weighs  15  pounds  and  was  made 
from  nonmagnetic  stainless  steel.  The  shaft  is  16  inches  long,  1-1/4  inches 
in  diameter  at  the  center  section  and  1/2  inch  in  diameter  at  the  ends.  The 
shaft  is  also  made  from  nonmagnetic  stainless  steel.  The  ratio  of  inertia  of 
the  rotating  assembly,  axial  vs.  transverse,  is  less  than  one  and  could  result 
in  stability  problems.  Preliminary  tests  at  low  speed  have  shown  that  this 
is  not  a problem  with  the  test  setup.  In  a real  system  this  ratio  will  be 
greater  than  one  and  thus  stable. 

The  stator  of  the  magnetic  bearing,  which  contains  the  Samarium  Cobalt 
(SmCO.)  magnets  and  the  coils,  is  mounted  to  the  fixed  base  structure.  The 
stator  is  made  from  high  silicon  low  chromium  steel  for  its  high  magnetir 
permeability.  The  magnetic  bearing  section  of  the  rotor  is  made  from  4130 
carbon  steel  . 

The  magnetic  gap  in  the  axdal  direction  is  0.015  inch  and  travel  is 
limited  in  both  the  axial  and  transverse  directions  by  the  touchdown  hall 
bearing  assemblies.  The  touchdown  bearings  engage  with  0.010  inch  travel 
and  protect  the  magnetic  bearings  from  rubbing  contact  within  themselves. 
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The  bill  bearing  assemblies  are  a duplex  pair  with  a twelve  pound  preLoad 
separated  by  a distance  of  10  inches.  The  inertia  of  the  rotating  element  of 
the  touchdown  bearing  assembly  was  minimized  so  that  its  acceleration  rate 
during  touchdown  at  high  speeds  will  be  very  high  and  thus  minimize  scuffing 
on  t lie  ba  1 Is. 

The  male  touchdown  stud  was  made  from  nonmagnetic  steel  and  the  rotating 
female  pieces  from  leaded  bronze.  The  bearings  were  high  speed,  duplex  turbine 
bearings  purchased  from  Barden  Corporation. 

The  geometry  of  the  touchdown  elements  are  bi-conical  at  45°  to  achieve 
touchdown  in  either  the  axial  or  transverse  directions  with  the  same  linear 
travel.  The  bi-conical  feature  of  the  touchdown  elements  limited  the  coning 
of  the  shaft  and  maintained  normal  alignment  at  the  extreme  axial  displacement. 
The  problem  of  fabrication  tolerances  of  the  touchdown  elements  was  controlled 
by  manufacturing  matched  sets  and  providing  for  shimming  at  assembly.  Assembly 
measurements  were  recorded  and  the  shims  ground  to  size  to  maintain  the  specified 
clearance  of  0.01  inch  within  1 .001  in. 

The  base  structure  consists  of  an  aluminum  bottom  plate  and  two  stanchions 
which  are  bolted  and  pinned  to  the  base.  The  stanchions  are  split  at  the  spin 
axis  to  facilitate  assembly  without  removing  the  stanchions  from  the  base.  The 
base  assembly  was  prematchined  and  then  assembled  and  final  machined  as  an 
integral  unit  to  maintain  alignment  of  the  axis.  The  base  structure  is  made 
from  6061 -Tb  aluminum. 

The  motor  generator  was  purchased  from  Sperry  Electro  Components,  Durham, 
North  Carolina.  The  unit  is  a brushless  DC  motor-generator  with  an  ironless 
armature  and  a 0.016  inch  axial  airgap.  The  motor-generator  consists  of 
three  components,  armature,  return  path  and  permanent  magnet  rotor.  The  arma- 
ture is  supported  to  the  base  structure  by  an  aluminum  bracket.  The  "return 
path"  and  permanent  rotor  are  supported  to  the  shaft  as  shown  in  Fig.  4-1. 

The  return  path  is  made  from  high  strength  stainless  steel.  The  permanent 
magnet  rotor  support  is  made  from  titanium. 

Based  on  dynamic  considerations  the  elements  of  the  assembly  were  designed 
for  high  stiffness  within  the  constraints  of  the  basic  configuration.  The 
desired  natural  frequency  of  the  Individual  elements  was  above  1 kHz,  a rather 
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ambitious  goal.  The  low  natural  frequency  of  the  shaft  (154  Hz)  was  a result 
of  placing  the  touchdown  bearings  inboard  of  the  magnetic  bearings.  It  was 
recognized  that  some  frequency  tuning  might  be  required  during  testing  to 
avoid  particularly  troublesome  frequency  bands. 

The  conventional  rotating  slender  shaft  configuration  is  not  the  optimum 
configuration  with  respect  to  stiffness  from  practical  considerations.  The 
fixed  shaft  configuration  offers  an  advantage  in  this  respect  and  should  be 
considered  in  advanced  designs. 

The  rotor  of  the  magnetic  bearing  assembly  must  react  a magnetic  force 
of  approximately  200  pounds.  The  deflection  of  the  rotor  due  to  this  force  is 
0.0006  inch.  Under  normal  operation  with  the  magnetic  bearing  centered,  the 
axial  forces  on  the  rotor  are  small.  The  axial  force  is  a function  of  the 
rotor-stator  clearance,  at  the  mid-position  the  force  is  small  and  increases  to 
a maximum  as  the  gap  is  decreased. 

The  stator  of  the  magnetic  bearing  assembly  is  inherently  stiff  with 
respect  to  both  the  inner  and  outer  pole  pieces  prior  to  assembly.  Fasteners 
were  not  permitted  in  the  assembly  of  this  unit  to  minimize  energy  losses. 

The  outer  pole  pieces  were  assembled  by  means  of  a shrink  fit  with  the  two 
structural  members. 

The  Samarium  Cobalt  (SmCO,-)  magnets  were  purchased  from  Spectraflex,  Inc., 
located  in  Watsonville,  California.  The  individual  magnets  were  premachined 
and  magnetized  prior  to  assembly  into  the  outer  pole  piece.  The  outer  pole 
pieces  were  fabricated  at  I.incoln  Laboratory  and  sent  to  Spectraflux  to  have 
the  magnets  installed.  After  bonding  the  magnets  to  the  pole  piece,  the 
inside  diameter  of  the  outer  magnets  was  ground  for  0.001  inch  diametrical 
clearance  with  the  inner  pole  piece.  The  combination  of  small  clearances, 
strong  magnetic  forces,  and  brittle  magnets  necessitated  the  use  of  sophisti- 
cated assembly  fixtures. 

Alignment  was  maintained  by  tolerance  control  of  machined  parts  and 
hand  fitting  at  assembly.  Spacers  were  provided  to  maintain  the  0.015  inch 
axial  gap  between  rotor  and  stator  of  the  magnetic  bearing  and  for  axial 
centering  of  the  touchdown  bearing  assembly.  Measurements  were  made  during 
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the  assembly  process  and  the  spacers  were  ground  to  achieve  the  specified 
clearances. 

A dynamic  analysis  of  the  energy  wheel  system  was  performed  to  determine 
resonant  frequencies  of  the  rotating  system.  A distributed  mass  was  modeled 
for  the  shaft  and  discrete  loads  for  the  flywheel  rotor  and  magnetic  bearing 
rotors. 

The  first  mode  frequency  of  the  rotating  element  was  calculated  to  be 
154  Hz.  The  first  mode  frequency  of  the  rotating  element  coupled  with  the 
magnetic  bearings  was  calculated  to  be  38  Hz.  Both  values  were  confirmed  by 
test . 

Pare  earth  magnets  exhibit  excellent  magnetic  properties  which  enable  an 
engineer  to  build  smaller  reliable  magnetic  parts.  The  common  permanent 
magnetic  materials,  alnlcos  and  ferrites,  are  insufficient  for  our  application 
because  of  size  and  weight  constraints,  and  reliability.  Samarium  Cobalt, 

SmCO^,  magnets  are  small  with  high  stability  and  provide  high  energy  which 
makes  them  adaptable  to  be  used  in  the  magnetic  bearing  In  the  energy  storage 
flywheel. 

All  magnets  Irrespective  of  the  manufacturing  techniques  have  fairly  well 
specified  remanence  and  induction  coercive  force;  however,  the  intrinsic  coercive 
force  and  the  energy  product  values  vary  over  a wide  range  as  a result  of  sinter- 
ing, heat  treatment  and  the  particle  size.  Magnetic  pieces  from  a single  lot 
cannot  be  accepted  because  several  test  samples  happen  to  meet  specifications. 
From  Table  1 it  can  be  seen  that  after  measuring  the  residual  inductance,  P.  , 
and  exposing  it  to  a 1()7.  demagnetization  knock  down  field,  H^,  487  failed  to 
meet  the  specification  requirements.  The  finished,  B^,  should  be  within  + 57 
of  the  lot  value  after  the  demagnetization.  Absolute  Br  could  not  be  readily 
measured  on  the  finished  product  because  the  shape  of  the  part  did  not  match 
the  induction  coil  utilized  in  the  measuring  procedure.  To  get  around  this, 
a standard  relative  to  the  experimental  setup  was  used  and  the  Br  after 
knockdown  was  compared  to  the  B^  before  for  each  piece. 

Properties  of  the  magnet  such  as  the  B^,  coercive  force,  Hc>  and  the 
intrinsic  coercive  force,  H|c»  are  merely  endpoints  in  the  intrinsic  or  the 
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inductive  demagnetization  curve  and  do  not  give  sufficient  information  about 
the  loop  shape.  Therefore  we  specified  that  a plot  of  the  second  quadrant  of 
the  magnetization  curve  be  supplied.  The  magnets  were  serialized,  permitting 
us  to  select  matched  sets  for  each  bearing. 

The  data  required  on  the  Samarium  Cobalt  magnets  that  we  considered  to  be 
most  useful  were  as  follows: 

«=  8,000  to  8,400  gauss 

H » 7,700  to  8,100  oersteds 

c 

H . > 18K  oersteds 

ci 

BH  > 16  * 10^  gauss  - oersteds 

and  full  second  quadrant  B-H  curve. 

The  high  silicon.  Si,  low  chromium,  Cr,  steel  used  in  the  stator  was 
supposed  to  give  a high  permeability,  y,  and  resistivity,  p.  Permeability 
expresses  the  dynamic  relationship  between  B and  H and  is  a function  B,  H,  p 
frequency  of  excitation,  etc.  Measurement  of  y and  p was  made  at  Lincoln 
Laboratory  and  they  fell  short  of  our  expectation.  However,  wet  chemical 
analysis  of  the  steel  gave  1.7%  Si  instead  of  about  37  but  an  emissl  n spectro- 
analysis  confirms  that  Cr  was  low  (<  0.5%). 

B.  Calculations 

The  target  design  of  the  brassboard  system  was  to  support  a 25-lb  rotor 
with  spin  axis  within  30°  of  vertical  in  1C,  and  allow  an  750  ft.  lb.  sec.  wheel 
to  be  rotated  at  approximately  0. 5°/sec.  cross  axis  rate.  This  was  derived  from 
system  studies  and  satellite  maneuvering  or  tipoff  requirements.  The  resulting 
design  impact  was  that  bearing  area  must  be  sufficient  to  supply  15  lbs.  of  lift 
(20%  safety  margin)  and  the  bearing  separation  adequate  for  the  cross-axis  torque. 

The  MTI  report  "A  Theoretical  and  Experimental  Investigation  of  the  Magnetic 
Fields  and  Forces  Arising  in  Magnetic  Suspension  Systems"  by  J.  Walowit  et  al., 
1975,  was  used  as  a basis  for  the  design. 

The  important  results  of  this  study  were  the  following: 


59 


1)  Optimum  tooth  width  T is  roughly  equal  to  the  pap  spacing  h 

2)  Optimum  tooth  separation  is  between  1 and  4 times  the  gap  spacing  h 

1)  Use  the  B field  in  the  middle  of  a tooth  B as  the  scaling  B field, 

o 

the  force  density  <*  1/2  h()  /Mn  38  the  scaling  force,  density,  the 

scaling  magnetic  flux  A » B A,  where  A is  the  tooth  area,  and  the 

o o t t 

MMF  **  BQh/Mo  where  h is  the  gap  spacing,  and  the  scaling  permeance 

is  A /0  - F 

o o o 

4)  With  the  above  scaling  values  and  optimum  dimensions,  for  square  teeth 

P/P  = 1.4,  F*/F  - 1.2,  F*/P*2  =0.6  F /P?  where  F*  is  the  force 
o o o o 

normal  to  the  tooth,  and  F( /F^  = 0.08  e/T  where  e is  the  t-ansverse 

displacement  and  F^  the  transverse  force  density  per  unit  total  bear- 
ing area.  Also,  F^  begins  to  flatten  out  for  e/T  > 0.5. 

Therefore,  T chose  to  allow  a reasonable  touchdown  bearing  clearance  of 
± (10  mils),  h = 15  mils,  T = 15  mils,  S = 50  mils.  I let  the  tooth  depth  be 
equal  to  the  spacing.  This  allowed  a working  bearing  range  of  i 7.5  mils  with- 
out saturating  F . Probable  value  of  the  bearing  force  density  would  at  e/T  = 

.5  then  be  0.04  F /2,  the  last  factor  of  two  due  to  the  cylindrical  geometry. 

The  above  choices  were  verified  by  comparison  with  available  designs  in  the 
literature.  Measured  F appears  to  be  10  to  20"'  higher. 

The  axial  forces  can  be  calculated  for  a single  model  of  our  bearing  with- 
out teeth  for  maximum  magnet  flux  A,  current  NI  in  the  coils,  pole  area  A, 

effective  total  gap  0^  / 2 per  pole),  and  displacement  x the  axial  force 

should  be 

_ 2 A NI  2 A2  x 

x 0 poA  G 

o o 

Therefore,  for  x = G /4  the  total  force  is  lust  equal  to  four  times  that 

° 2 

calculated  for  one  unbalanced  pole  of  po  /8  A,  and  the  force  per  rotor  equal 
to  the  unbalanced  equilibrium  force  on  one  rotor. 
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So,  to  estimate  axial  stiffness,  at  B()  = .8')  T the  average  flux  will  be 
! . 7B((  times  the  tooth  area  or  .33  T per  unit  pole  area.  This  gives  an  unbal- 
anced pole  force  density  of  (.33/. 83)^  F = .156F  In  equilibrium.  Therefore, 

o o 

at  7.3  mils  displacement  we  will  have  F = .156F  and  R„  » .02F  for  a stiff- 

A o K o 

ness  ratio  of  7.8,  which  is  In  good  agreement  with  the  literature  values  of  6 
to  8 . 

For  most  kinds  of  Iron,  tooth  saturation  occurs  at  B = 2.0  T and  one  must 
allow  perhaps  a factor  of  two  due  to  field  concentrations  near  corners.  The 
desired  field  F was  set  to  be  .0833  T.  This  gave  F = 2.76  x 10^  n/m2  and 

3 o 2 o 

F(  = 5.5  x 10  n/m  per  area  of  pole  In  the  transverse  directions. 

The  total  30  Ibf  = 133  N is  supported  effectively  at  1/2  G at  two  ends, 

giving  33.25  N per  bearing,  16.6  per  rotor,  or  8.3  per  pole.  Therefore,  the 

2 2 
pole  area  must  be  15.1  cm  , or  a tooth  area  3/13  of  this  = 3.48  cm  of  tooth 

area,  including  the  20%  safety  factor  of  force  required. 

Based  on  a shaft  diameter  of  .75  inch,  an  inner  radius  for  the  bearing- 

pole  pieces  of  1 cm  was  required.  Then,  by  adding  up  tooth  areas  with  15  mil 

thickness  and  50  mil  spaces  between,  the  radius  of  the  inner  teeth  could  be 

terminated  at  9 teeth  with  2.285  cm  radius  and  3.485  cm2  accumulated  area. 

heaving  just  over  1 cm  for  colls,  the  inner  ring  of  the  outer  poles  had  a 

radius  of  3.30  cm.  Adding  rings  and  spaces  as  before  give  four  teeth  with 

2 

3.83  cm  radius  = 1.5  inches,  with  another  3.38  cm  accumulated  area  (close 
enough) . 

Now  to  determine  the  required  spacing  between  bearings.  At  15  lbs  with 

motion  of  7.5  mils,  the  radial  stiffness  is  K = 2000  Ibf/in.  and  for  our  stiff- 

r 

ness  ratio,  the  axial  stiffness  is  = 16,150  lbf/in.  If  all  rings  were 
located  at  a single  radius  D and  separation  L from  the  center  of  mass,  the 
torsional  stiffness  would  be 


kt  = krl2  (i  - kad2/krl2) 

Suppose  (worst  case)  I)  - 1.5"  and  L » 6.0".  Then, 
in-lb/rad  or  3028  ft-lb/rad. 


.252  K^L.2 


3.64  x 10 


4 


61 


With  angular  deviation  of  7.5  mils/6  in.  c .00125  rad.  This  allows  3.8  ft- 

_3 

lb.  or  a cross-axis  rate  of  5 x 10  rad/sec  ■ .28  deg/sec.  However,  things 

are  not  this  bad  since  most  rings  are  located  closer  to  the  spin  axis,  the 

average  radius  (center  of  the  coils)  being  at  1.1  in.  so  that  a better  estimate 

2 

of  torsional  stiffness  is  K *=  .37  K I,  and  the  resulting  cross-axis  rate  capa- 

T K 

bility  is  .413  deg/sec.  This  is  close  enough  to  specified  values  for  a first 
cut  design,  since  the  relative  stiffnesses  are  only  approximately  known.  If 
the  stiffness  ratio  is  actually  at  the  low  end  (6.0)  this  gives  the  required 
cross-axis  torque. 

To  achieve  the  field  of  Bq  - .833  T in  the  center  of  a tooth,  assuming 
average  B/B^  is  perhaps  0.6,  the  required  A/A  » 1.75  so  the  average  B in  the 
iron  behind  the  teeth  is  1.75  (3/13)  Bq  * 0.33  T.  The  iron  can  be  used  to  con- 
centrate flux  at  the  permanent  magnet  to  be  0.67  T,  about  80%  of  the  residual 
flux  level  of  SaCo  magnets. 

TJhe  required  flux  from  the  magnet  is  double  the  1.75  B times  the  tooth 

2 0 
area  of  3.48  cm  , or  10.1  webers.  At  the  flux  level  of  .67  T this  requires  a 

2 

magnet  cross  section  of  15.1  cm  , at  an  average  radius  of  3 cm.  This  then 
requires  an  axial  length  of  0.8  cm.  Allowing  20%  leakage  flux  gives  a cross 
section  requirement  for  a magnet  ring  at  1.2  inches  radius  with  length  of 
0.4  inch. 

The  expected  MMF  requirement  is  2(0.8)  Bq  h/mQ  or  400  amps.  At  70%  of  Br 
a SaCo  magnet  can  be  expected  to  give  an  MMF  of  about  20%  of  its  coercive  force, 
or  1600  oersteds  of  1.27  x 10"*  amps/meter.  This  then  requires  a magnet  thick- 
ness of  3.17  mm  - 125  mils.  To  give  a safety  factor  and  allow  for  iron  losses, 
but  primarily  to  increase  magnet  thickness  for  physical  integrity  with  a better 
form  factor,  this  was  Increased  to  150  mils. 

To  keep  drag  low,  magnet  sections  are  to  be  matched  within  5%  or  better 
and  taken  from  a single  lot.  Also,  bearing  tooth  dimensions  are  to  be  tightly 
controlled  to  determine  how  low  the  drag  can  be  pushed. 


62 


C . Loop  Design 


In  the  axial  direction,  the  stiffness  of  the  bearing  will  be  about  8x 

3 A 

the  passive  direction  stiffness  of  2 x 10  lbf/in.,  or  1.6  x 10  lbf/in.  = 

2.8  x 10^  n/m.  The  rotor  weighs  about  25  lbtn  - 11.34  kg.  Therefore,  the  axial 
resonant  frequency  is  497  rad/sec  • 79  Hz. 

The  loop  compensation  was  intended  to  create  a crossover  frequency  lOx 
higher  than  this,  with  a midband  gain  of  30  and  a phase  margin  of  45°.  To  do 
so,  a zero  was  placed  at  3x  the  axial  resonant  frequency  and  a pole  at  30x  that 
frequency.  This  should  give  the  45°  desired  phase  lead  at  a crossover  of  800  Hz. 
All  other  poles  must  be  located  well  beyond  2.4  kHz,  preferably  beyond  8 kilz. 

Therefore,  a sensor  with  10  kHz  bandwidth  has  been  chosen.  Effects  of 
structural  resonances  and  phase  lag  due  to  eddy  currents  in  the  bearing 
materials  must  be  measured  and  compensated  for  during  system  testing. 

During  integration,  it  was  found  that  the  many  vibration  modes  of  the  system 
required  lowering  the  crossover  frequency  to  150  Hz. 

D.  Magnetic  Bearing  Integration  and  Test 

The  assembly  of  magnetic  bearing  parts  went  relatively  smoothly.  Press- 
fitting of  stator  parts  required  good  positive  control  of  the  planar  positioning 
of  parts;  otherwise  they  tended  to  tip  over  and  bind.  One  stator  was  success- 
fully disassembled  and  rebuilt. 

Stators  were  measured  for  magnetic  field  variation  with  rotation.  A mod- 
ification of  the  fixturing  for  stator  assembly  and  a magnetoresistor  were  used. 
The  results  showed  that  the  open-circuit  stator  field  of  1000  gauss  had  varia- 
tions on  S/N  1 of  + 20  gauss.  These  variations  should  be  halved  again  at  least 
by  the  addition  of  rotors  to  close  the  circuit.  Therefore,  variations  in  B of 
0.2%  can  be  expected  with  rotation. 

Mechanical  tolerances  on  bearing  tooth  runout  and  height  were  measured  to 
be  typically  + 0.1  mil,  with  occasional  peaks  of  ± 0.2  mils.  These  represent 
about  0.3%  variations  in  the  effective  gap  total  of  30  mils. 

Therefore,  the  sum  of  all  B-field  variations  1b  expected  to  be  within  0.5%. 
This  will  lead  to  drag  due  to  eddy  currents  in  the  bearing  material  with  power 


63 


1 2 

density  of  — (oiABL)  where  L is  the  length  scaling  the  eddy  currents.  Over  a 

12  2 3 

volume  = AL  the  resulting  power  loss  is  on  the  order  of  P = -r—  to  AB  L A watts. 

-5  -7  -A  -3  2 

For  measured  p = 10  fi-cm  = 10  12-m,  L = 3.8  x 10  m,  A = 4 x 10  m, 

A3  = 3.5  x 10  ^T,  w = 5 x lo\  P = 3 x 10  ^ watts.  Ridiculous!  But  if 

_2 

L = 3.8  x 10  m,  then  P = 300  watts,  so  the  eddy  current  scaling  length  L is 
critical . 

Assembly  of  the  bearings  onto  the  shaft  went  as  planned.  Forces  near  250 
pounds  were  encountered  when  positioning  bearing  pole  pieces  and  stators  onto 
the  shaft,  but  the  use  of  jacking  screws  and  holding  fixtures  was  successful. 

Use  of  aluminum  or  stainless  in  fixtures  caused  some  galling  problems,  which 

were  solved  by  a little  lubricant.  A better  answer  would  be  to  use  brass  for 
tight  fitting  fixtures. 

Applying  swept  frequency  currents  to  the  motors  allowed  measurement  of 
resonances  in  the  structure.  Also,  the  current  needed  for  moving  the  rotor  from 
one  end  to  the  other  was  measured  to  be  2 amps  on  both  motors,  or  4 amps  on  one 
motor.  This  showed  some  increase  (to  near  6 amps  on  one  end)  as  time  progressed, 
perhaps  due  to  extended  touchdown  gaps.  The  linearity  was  gratifying. 

Resonances  were  found  at  several  frequencies  in  the  structure.  One  problem 
was  the  original  position  sensor  and  its  mounting  bracket.  It  had  internal 
nodes  in  the  micrometer  adjustment  as  well  as  a cantilevered  geometry.  A new 
sensor  with  no  micrometer  and  a bracket  with  no  cantilever  was  designed  and 
solved  this  problem. 

The  flywheel  rotor  turned  out  to  have  a very  high  Q resonance  at  about 
2.1  kHz.  This  was  clearly  audible  as  a ringing  if  the  rotor  were  struck.  Use 
of  fiberglass  spacers  between  rotor  and  shaft  did  not  help.  Eventually,  a notch 
filter  with  Q = 5 was  built  to  eliminate  the  problem.  This  worked.  A non- 
metallic  rotor  would  eliminate  the  problem. 

Several  other  resonances  were  not  easy  to  locate  but  did  not  cause  large 
phase  shifts  in  the  bearing  loop.  A transverse  resonance  due  to  shaft  bending 
at  154  Hz  was  numerically  predicted  and  verified.  Longitudinal  resonance  at 
400  Hz  and  at  600  Hz  is  believed  to  be  bending  of  the  yoke  arms.  A stiffening 
yoke  member  may  be  added  later. 
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The  bearing  was  floated  initially  with  the  new  sensor  and  bracket  but  no 
notch  filter. 

It  had  no  phase  margin  and  a small  (±1  mil)  residual  oscillation  at  its 
crossover  frequency  of  150  Hz.  Addition  of  the  notch  filter  allowed  lead 
compensation  which  gave  over  30°  of  phase  margin,  a dc  loop  gain  of  5,  and 
crossover  frequency  of  150  Hz.  The  bearing  is  sensitive  to  resonances  in  its 
mounting  structure.  Mounting  on  a table  which  is  not  flat  causes  oscillations 
at  about  400  Hz,  believed  to  be  due  to  the  rocking  of  the  assembly  causing 
phase  shifts  and  interacting  with  the  yoke  bending  modes. 

In  use  in  a satellite,  such  a bearing  would  have  to  be  designed  into  the 
structural  dynamics,  and  the  dynamics  of  the  structure  also  designed  into  the 
bearing  loop. 

The  bearing  motor  inductance  was  close  to  the  15  mH  predicted  value  at  low 
frequencies  and  as  predicted  fell  off  as  almost  l//s  at  frequencies  above  100  Hz 
due  to  eddy  currents.  To  avoid  large  circulating  currents  of  50  kHz  due  to 
current  driven  switching,  a 2 mH  pot  core  inductor  was  added  in  series. 

The  floated  bearing  showed  stiffnesses  of  about  10  lbf/mil  along  the  spin 
axis,  about  1.2  lbf/mil  transverse,  and  3.26  foot-lbf  per  milliradian  transverse 
rotational  stiffness.  This  was  in  good  agreement  for  axial  to  radial  stiffness 
but  low  by  about  a factor  of  two  in  total  stiffness.  The  cause  is  presently 
unknown.  Increased  stiffnesses  can  probably  be  obtained  by  reducing  the  magnetic 
gap. 

One  possible  explanation  of  the  low  stiffness  values  could  be  larger  than 
predicted  reluctance  in  the  bearing  pole  piece  iron.  This  could  result  in 
increased  flux  leakage  or  in  the  flux  clustering  near  the  magnet  and  not  giving 
much  centering  at  the  inner  and  outer  fringe  rings.  Such  an  effect  could  be 
avoided  by  designing  for  larger  magnetic  gap  and  larger  clearances  (more  magnet 
length),  or  by  using  better  iron  or  lower  field  strengths  (more  bearing  area). 

Due  to  the  small  dimensions  of  the  gap,  it  is  difficult  to  measure  the  radial 
field  distribution. 

Resonant  frequencies  were  looked  for  to  verify  stiffness  measurements  and 
check  damping.  The  transverse  translation  resonance  at  22  Hz  checks  out,  and 
has  a damping  time  constant  of  about  2 seconds.  The  rotational  damping  is  too 
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large  to  measure;  it  may  be  critically  damped.  A Bode  plot  of  position  vs. 

Input  current  shows  a very  linear  f orce/current  relationship  and  a pair  of 
poles  at  80  Hz  as  expected. 

Position  sensor  outputs,  using  Lion  precision  capacitance  non-contact 
sensors,  are  used  for  position  and  velocity  sense  lines.  These  are  very 
troublesome  due  to  the  capacitance  sensitivity  of  all  1 MHz  tuned  cables  and 
non-interchangeable  probe  types  causing  anomalously  high  gain  (almost  double) 

In  the  sensor  outputs.  These  are  supposed  to  be  approximately  1 volt/mil  and 
1 vol t /m il/mill lsecond , respectively.  The  small  signal  electronics  transfer 
function,  without  the  2.1  kHz  center  Q:5  notch  filter,  is  then 

^(amps)  (2  ms  S+l) 

Z(mils)  “ (.039  ms  S+l)  (.047  ms  S+l) 

(1  ms  S+l)  (.33  ms  S+l)  (.07  ms  S+l)  (.02  ms  S+l) 

(.5  ms  S+l)  (.14  ms  S+l)  (.04  ms  S+l)  (.22  ms  S+l) 

The  alternating  poles  and  zeros  in  the  right  side  of  this  expression  are  an 
attempt  to  compensate  for  the  1/ Js  eddy  current  effect  with  a 45°  lead  network. 
The  zero  at  2 ms  is  created  by  adding  the  position  and  velocity  outputs  of  the 
sensor  together  with  appropriate  gains.  The  poles  at  .039  ms,  .047  ms,  and 
.22  ms  are  low-pass  filters  to  eliminate  the  effects  of  high-frequency  reso- 
nances In  the  system. 

An  additional  feedback  loop  using  an  integrator  with  clipped  output  is  used 
to  eliminate  steady-state  current  by  automatically  seeking  the  null  force  point. 
This  "VZP"  circuit,  first  used  in  this  application  by  S.  Lyman  of  Cambion  Corp., 
adds  a multiplier  to  the  transfer  function 

100  ms  S 
100  ms  S-l 

With  this  circuit  in  operation,  or  without  it  and  with  the  offset  position 
trimmed  to  zero,  the  power  drawn  by  the  fully  redundant  bearing  system  is  just 
under  2.0  watts.  This  is  independent  of  the  rotor  weight  up  to  the  touchdown 
point:  over  150  pounds  in  the  vertical  position.  This  input  power  can  probably 
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be  reduced  by  a factor  of  two  In  future  systems  if  it  matters  by  improving  tbe 
sensor  electronics  and  reducing  the  number  of  op  amps.  Redundancy  has  been 
demonstrated  by  manually  switching  off  one  motor  and  turning  on  the  other,  one 
sensor  per  motor. 

To  enhance  the  liftoff  (large  signal)  stability  of  the  bearing,  a nonlinear 
lead  network  is  added.  This  is  to  attempt  to  compensate  for  the  I.dl/dt  limita- 
tions of  the  power  supply  (+  15  volts).  Through  Xener  diodes,  the  network 
introduces  a 4.7  ms  zero  above  2.4  amps  and  increases  this  to  a 22  ms  zero  above 
5 amps,  with  an  accompanying  .044  ms  pole. 

The  current  drive  circuit,  presently  a half  bridge  switched  at  50  kHz,  is 
quite  efficient  but  should  be  replaced  by  a full  bridge.  This  would  allow 
larger  Ldl/dt,  would  simplify  power  conditioning  by  using  unregulated  bus,  and 
would  avoid  one  supply  boosting  the  other  through  the  current  driver  output 
inductance . 

The  motor  characteristics  of  the  bearing  are  50  Ibf/amp  and  5 lbf/mil  per 
motor  giving  a measured  total  open  loop  dc  gain  of  5 for  one  active  motor  and 
one  passive  one.  If,  as  appears  possible,  the  sensor  gain  is  anomalously  high 
X2 , the  motor  gain  may  be  lower  by  a factor  of  two. 

Magnetic  gaps  have  been  measured  in  floated  and  non-floated  conditions  to 
be  within  5%  of  nominal  values:  gaps  on  both  sides  of  the  stator  sum  to  30  mils. 
Total  rotor  motion  between  touchdown  stops  is  22.5  mils.  Drag  is  dominated  by 
air  flows  and  spin-down  time  constant  is  measured  to  be  one  to  two  hours  in  air 
at  ambient  pressure.  Ball  bearing  drag  by  comparison  gives  spin-down  time 
constants  of  less  than  one  minute  under  150  lb.  estimated  preload  on  the  touch- 
down bearings. 

The  actual  hardware  under  test  is  shown  in  Fig.  4-2.  Construction 
details  are  shown  for  the  stator  in  Fig.  4-3,  4-4,  and  4-5,  the  test  flywheel 
and  shaft  in  Fig.  4-6,  the  touchdown  bearings  in  4-7,  bearing  housing  and  shims 
in  4-8,  assembly  fixtures  in  4-9  and  the  electronics  control  panel  in  4-10. 

Several  very  interesting  observations  were  made  while  measuring  bearing 
drag.  First,  the  drag  is  dominated  by  aerodynamic  forces  down  to  quite  low 
pressures,  on  the  order  of  10  microns.  This  Implies  that  a spacecraft  device 
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Fig.  4-3.  Bearing  stator  assembly  and  test  bed 


Fig.  U-U.  Construction  of  bearing  stator. 
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Fig.  4-7.  Touchdown  bearings. 
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must  be  well  vented  to  take  full  advantage  of  its  environment,  and  ground  test- 
ing should  be  done  with  a cryopump.  Second,  the  drag  has  its  lowest  value 
below  the  lowest  critical  speed,  reaches  a sharp  maximum  at  that  speed,  and 

settles  down  to  a higher  drag  at  supercritical  speeds.  For  our  test  system, 

2 

rotating  M = 25  lbs.  and  I = 0.586  in.lbf-sec  . Besides  the  resonance,  the 

drag  torque  appeared  to  be  linearily  proportional  to  speed  over  a range  of 

50  RPM  to  1.5  KRPM,  with  critical  speed  at  about  1.2KRPM.  Below  critical, 

if  T = Bw  and  spindown  time  constant  x = I/B  = 30  hours  and  so  B = 5.4  x 10  ^ 

-5  -3 

in-lbf-sec  = 8.7  x 10  in-oz-sec.  or  a drag  of  9.1  x 10  in-oz.  per  1KRPM. 

-4 

Above  critical  speed,  drag  was  approximately  22  hours,  or  B = 1.2  x 10 
in-oz-sec.  for  a drag  of  .012  in-oz.  per  1KRPM.  This  was  measured  with  a rotor 
which  had  not  been  balanced  at  all.  Therefore  it  should  not  be  any  surprise 
that  the  drag  is  higher,  due  to  the  rotation  of  the  shaft  geometric  center 
about  the  principal  axis  of  the  rotor.  In  fact,  this  indicates  that  poorly 
balanced  nonmetallic  rotors  will  not  have  a problem  due  to  drag. 

Drag  versus  speed  and  pressure  are  plotted  on  Figs.  4-11  and  4-12.  The 
30-hour  drag  rate  corresponds  to  power  inputs  of  .67  watt  at  10,000  RPM  or 
19  watts  at  50,000  RPM. 

E.  Finite  Element  and  Field  Plotting  Capabilities 

To  allow  exploration  of  magnetic  bearing  tooth  geometries  and  to  verify 
design  calculations,  a magnetic  field  finite  element  analysis  program  has  been 
developed.  This  was  used  to  verify  reluctance  calculations  for  the  LES  magnetic 
bearing.  Further  work  using  this  program  is  expected  to  give  more  insight  into 
construction  of  many  types  of  electromagnetic  devices. 

MIT/LL  Electromagnetic  Field  Computer  Analysis 

Lincoln  Laboratory  has  added  electromagnetic  field  calculation  features 
to  the  well-known  general-purpose  structural  analysis  system  ICES-STRUDL.  When 
operated  on  the  Laboratory's  IBM  370/168  computer,  the  ICES  executive  provides 
the  data  management  and  program  management  facilities  needed  to  minimize  stor- 
age and  computation  requirements,  while  the  STRUDL  subsystem  provides  the 
finite  element  and  equation  solving  facilities  needed  for  efficient  analysis. 
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Families  of  isoparametric  finite  elements  are  available  for  scalar  poten- 
tial models  (electric  or  magnetic)  and  for  magnetic  vector  potential  models. 

These  families  presently  consist  of  6 different  types:  linear,  quadratic, 
and  cubic  planar  quadrilaterial  elements;  linear  and  quadratic  3-dimensional 
brick  elements;  and  a linear  3-dimensional  triangular  prism  element.  Nonlinear 
material  properties  for  permeability  and  permittivity  will  soon  be  available. 

Output  quantities  available  are  scalar  or  vector  potentials,  flux  densities, 
field  intensities,  energy  densities,  and  electromagnetic  forces.  These  results 
may  then  be  graphically  interpreted  using  the  general-purpose  structural  analy- 
sis plotting  system  ICES-STRUPLOT.  The  STRUPLOT  subsystem  operates  in  the 
interactive  time-sharing  environment  of  CMS-ICES,  using  free-format  CDL 
(command  definition  language).  Commands  are  selected  dynamically  and  are  based 
on  visual  observation  of  successive  plots  produced. 

Enclosed  are  four  typical  ICES-STRUPLOT  plots  from  a magnetic  bearing  tooth 
case  study  using  quadratic  eight-noded  planar  elements.  The  plots  verify  the 
geometrical  accuracy  of  the  finite  element  modeling  (Fig.  4-13);  illustrate  the 
contour  plotting  feature  for  potentials,  field  intensities,  and  flux  densities 
(Figs.  4-14,  4-15,  4-16);  and  illustrate  nodal  averaging  for  smoothing  results 
(Fig,  4-16).  Note  that  the  plots  are  a selected  subset  of  the  total  model  and 
also  that  specific  components  of  the  field  intensity  and  flux  density  vectors 
can  be  selected. 

These  described  electromagnetic  finite  element  and  field  plotting  capa- 
bilities are  an  extension  of  well-known  and  popular  general-purpose  structural 
analysis  systems,  and  future  maintenance  and/or  further  enhancement  of  them 
will  require  very  little  extra  effort. 
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V . MOTOR  GENERATOR 


A.  Design  Goals  and  Tradeoffs 

An  efficient  flywheel  system  needs  a motor  generator  with  very  low 
coasting  losses  and  efficient  generating  and  motoring  modes,  in  that  order  of 
importance.  High  speeds  allow  a smaller  unit  and  smaller  electronics  with 
higher  quality  output,  but  make  the  efficiency  problem  worse.  By  using  uncon- 
ventional design  approaches,  the  LL  motor  generator  program  attempted  to 
resolve  this  difficulty.  The  attempt  appears  to  be  successful  based  on  avail- 
able diagnostics.  High  speeds  require  no  slip  rings  or  moving  contacts.  It 
is  difficult  but  not  impossible  to  keep  field  variations  small  with  a homo- 
polar  or  Lundell  type  induction  machine,  so  that  a permanent  magnet  rotor 
looked  attractive.  Active  solid  state  commutation  would  allow  output  regu- 
lation, so  the  lack  of  control  of  field  strength  would  not  be  a problem. 

Also,  the  lack  of  field  exitation  would  avoid  associated  losses.  However,  the 
field  could  not  be  turned  off  during  coasting  to  avoid  drag. 

The  availability  of  Samarium  Cobalt  permanent  magnets  prompted  to  choice 
of  a permanent  magnet  rotor.  They  offer  a smaller,  lighter  rotor  than  would 
otherwise  be  possible,  and  allow  large  magnetic  gaps  to  minimize  armature 
losses . 

To  avoid  high  drag  from  the  rotating  field,  an  ironless  stator  was  chosen. 

This  is  compatible  with  the  high  coercive  force  of  the  rotor  magnets.  To  avoid 

eddy  current  losses  in  the  stator,  the  conductors  were  made  of  Litz  wire  (very 
small,  individually  insulated  filaments)  and  no  other  conducting  material  was 
used  in  the  rotor.  All  iron  in  the  magnetic  path  rotates  with  the  magnets. 

The  choice  of  numbers  of  poles  could  be  made  two  ways.  If  the  stator  cur- 
rents are  to  be  sinusoidal,  produced  by  a chopped  DC  waveform,  then  the  number 

of  poles  must  be  minimized  to  keep  the  chopping  frequency  low.  At  50  KRPM,  a 
two-pole  machine  with  ten  chopping  cycles  per  rotation  requires  over  8 kHz 
chopping  rate.  A two-pole  machine  would  require  long  wire  runs  between  poles. 

However,  using  the  concept  of  square-wave  drive,  the  chopping  frequency 
equals  the  electromagnetic  angular  rate,  so  that  ten  times  the  number  of  poles 
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can  be  accommodated  at  the  same  switching  speed.  It  is  also  desirable  to  use 
a 3 phase  drive  to  eliminate  third  harmonic  effects  in  the  square-wave  drive. 

A 12-pole  design  allows  this,  with  electrical  chopping  frequency  of  Ax  rota- 
tional rate,  which  at  50  KRPM  gives  about  3.3  kHz.  Allowing  10%  for  rather 
thin  pulses  for  low-power  motoring,  this  gives  pulse  widths  of  15  ysec,  com- 
fortably longer  than  the  1 ysec  switching  times  typical  of  appropriately 
sized  transistors. 

A scheme  for  laying  two  rows  of  conductors  with  all  three  phases  and  no 
empty  space  was  worked  out,  allowing  minimum  gap  in  the  motor  and  so  minimum 
magnet  size.  Stresses  in  the  magnets  appeared  critical  and  a system  of  support 
using  titanium  for  support  and  an  iron  pole  piece  was  designed.  An  iron  return 
path  rotor  was  used. 

The  resulting  structure  was  therefore  a SaCo  and  steel  and  titanium  rotor, 
an  iron  return  path,  and  an  ironless  stator  with  two  layers  of  Litz  wire. 

Twelve  poles  and  three  phases  were  selected. 

Various  electrical  voltage  level  schemes  and  switching  circuits  were  con- 
sidered. The  driving  consideration  was  that  generating  efficiency  be  high,  so 
the  open  circuit  voltage  was  chosen  to  be  higher  than  the  required  28  VDC  bus 
over  the  generating  range.  An  upconverter  to  a variable  DC  bus  is  used  for 
motoring,  and  a 3-phase  bridge  of  diodes  and  transistors  (6  each)  used  to 
interface  the  M/G  leads  to  a varying  DC  bus.  The  variable  DC  is  then  down- 
converted  (efficiently)  to  the  28  VDC  satellite  bus. 

At  this  point  in  the  design,  it  was  decided  to  contract  out  the  final 
design  and  characterization  of  the  motor-generator  to  an  experienced  manu- 
facturer. A design  specification  for  the  motor-generator  and  interface  draw- 
ing were  drafted  and  released.  Goals  of  5 lbs  and  85%  generating  efficiency 
were  set  for  a fully  redundant  250  watt  unit. 

B.  Design  Realization 

The  contract  was  awarded  to  Sperry  Electro  Components  after  competi- 
tive bidding.  Sperry  reviewed  the  design  tradeoffs  in  the  areas  of  chopping 
frequency,  conversion  circuitry,  and  hardware  implementation.  Their  findings 
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agreed  with  our  own  that  a three-phase,  square-wave  drive  should  be  used,  that 
12  poles  was  near  optimum,  that  the  bridge  and  motoring  upconverter  and  genera- 
ting downconverter  g ve  best  system  efficiency,  and  that  weight  and  efficiency 
goals  were  close  to  achievable  values.  The  titanium  magnet  support  structure 
was  modified  but  conceptually  adopted.  The  use  of  Litz  wire  was  adopted  too, 
but  more  to  improve  field  penetration  into  the  windings  than  to  decrease  drag. 
The  coils  were  wrapped  on  three  separate  layers  for  ease  of  manufacture,  with 
fiberglass  stiffness  between  coils.  The  connection  scheme  was  modified  from 
delta  to  wye  to  allow  for  nonuniformity  of  fields  without  excessive  losses. 
Appendix  A is  the  Sperry  report. 

Sperry  built  a brassboard  of  the  motor  generator.  It  is  like  a final  model 
except  for  the  use  of  readily  available  steel  instead  of  high-strength  materials 
as  it  is  unable  to  run  at  full  50  KRPM.  This  brassboard  can  be  mounted  to  the 
LES  magnetic  bearing  test  bed  for  further  experiments. 

Results  of  the  design  show  that  even  higher  efficiencies  than  expected  can 
be  achieved  with  this  unit,  with  average  constant-power  generating  efficiency 
of  the  generator  of  96%  and  electronics  92.5%,  giving  system  overall  generating 
efficiency  of  88.5%  estimated.  Motoring  is  also  good,  with  machine  94.5%  and 
electronics  89.5%  for  overall  motoring  efficiency  84%.  This  gives  total  input- 
output  efficiency  of  74.3%,  a very  respectable  number. 
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VI.  CONCLUSIONS 


Momentum  wheels  made  with  magnetic  bearings  are  far  superior  to  those  with 
ball  bearings  for  spacecraft  use.  They  allow  higher  speeds  and  therefore  are 
lighter.  They  are  more  reliable  over  long  periods  of  time.  They  require  very 
little  operating  power  and  can  be  made  fully  redundant.  Their  design  should 
be  radially-passive,  axially-active  for  applications  requiring  low  cross-axis 
rates  due  to  satellite  maneuvering,  and  5-axis  active  for  applications  requiring 
higher  (over  1 or  2 degree  per  second)  cross  ayis  rates  or  an  electronically 
gimballed  wheel. 

The  most  troublesome  aspect  of  a momentum  wheel  with  magnetic  bearings  is 
the  potential  interaction  of  the  active  bearing  loops  with  the  satellite  struc- 
ture. This  can  be  resolved  by  design  analysis  but  must  not  be  neglected. 

Energy  storage  using  flywheels  is  only  slightly  lighter  than  NiH^  batteries 
and  is  more  complex.  However,  such  energy  storage  looks  attractive  if  atti- 
tude control  is  combined  with  it.  It  also  looks  attractive  if  power  require- 
ments dictate  extensive  conditioning  of  the  battery  output  voltage,  through  a 
regulator  or  inverter. 


APPENDIX  A 


Lincoln  Motor -Generator  Design  (Contract  with  Sperry  Rand  Corporation) 

The  Lincoln  Laboratory  flywheel  system  motor-generator  design  is  compatible 
with  system  design  layout  Fig.  4-1.  It  is  a fully  redundant  250  watt  DC  nrush- 
less  unit  for  operation  on  a 28  VDC  satellite  bus.  Its  samarium-cobalt  perman- 
ent magnet  rotor  is  capable  of  50,000  RPM  operation  with  full  generation  capa- 
bility down  to  20  KRPM.  Extra  windings  are  used  for  induced  EMF  commutation 
with  a special  start-up  circuit  needed  for  cold  starting.  It  is  a 3-phase, 

12  pole,  ironless  armature  design.  The  three  phases  are  separately  wound  on 
individual  non-conducting  armature  stiffeners  for  ease  of  assembly,  using  Litz 
wire  to  minimize  eddy  current  losses.  The  permanent  magnets  and  iron  in  the 
magnetic  circuit  all  rotate  together  to  eliminate  eddy  currents  in  the  iron. 

A titanium  containment  structure  houses  the  magnet9  and  iron  on  the  rotor  to 
support  the  magnet  material  at  high  speeds.  Magnets  are  purely  in  compression 
by  design.  Estimated  efficiency  as  a generator  is  over  90%,  as  a motor  over 
85%,  including  electronics.  Performance  analysis  and  brassboard  hardware  were 
delivered  on  1 October  1977. 

This  design  should  scale  up  easily  to  the  6-10  kW  level  for  larger  system 
applications.  Higher  scaling  levels  might  impact  the  power  switching  electronic 
interface. 

Sperry's  report  on  the  Motor  Generator  contract  follows. 
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A . 1 INTRODUCTION  AND  PURPOSES 


This  report  is  written  to  summarize  the  work  performed  on  Massachusetts 
Institute  of  Technology  Lincoln  Laboratory  Purchase  Order  F19628-76-C-002  PO 
CX-1340,  Flywheel  Motor-Generator  Design  as  covered  by  MIT  Specification  dated 
8 February  1977  and  RFI  No.  30098  dated  21  December  1976  and  subsequent  modifi- 
cations as  discussed  in  the  design  coordination  meetings. 

This  order  covered  the  design  of  an  Axial  Airgap,  Ironless  Armature, 
Brushless  DC  Motor-Generator  and  the  development  techniques  for  its  fabrication. 
The  design  effort  involved  trading-off  motor  generator  performance  with  the 
driver  type  and  then  optimizing  the  generating  performance  with  secondary  goals 
of  minimum  weight,  size,  and  high  motoring  efficiency. 

An  Engineering  sample  was  fabricated  to  verify  manufacturing  techniques 
and  as  a check  on  the  magnetic  circuit  design. 

A. 2 REQUIREMENTS 

The  primary  goal  of  this  design  effort  was  to  obtain  maximum  efficiency  in 
the  generator  mode.  Additional  requirements  and  design  goals  were  used  to  con- 
strain the  design.  The  salient  requirements  for  the  Motor-Generator  are  pre- 
sented in  Table  A-l. 

A. 3 SUMMARY  OF  RESULTS 

A tradeoff  analysis  was  made  of  several  possible  commutation/regulation 
systems  resulting  in  selection  of  a three-phase,  discretely  commutated  motor 
generator  using  an  external  regulator  and  operating  at  relatively  high  voltage, 
i.e.,  always  higher  than  the  bus  voltage. 


A twelve-pole  motor/generator  was  designed  which  exceeded  all  performance 
requirements.  This  design  has  the  following  features: 


Fully  Redundant  Windings. 

Samarium  Cobalt  Magnets. 

Integral  Plckoff  Winding  for  Commutation. 

Primary  Losses  Resulting  In  96%  Generator  Efficiency. 

Lltz  wire  for  Coils  to  Minimize  Losses. 

High  Strength  Titanium  Rotor  Structure. 

Figure  A-l  shows  the  performance  and  outline  drawing  for  this  unit. 

A. 4 SELECTION  OF  MOTOR/GENERATOR  TYPE  AND  DRIVE  INTERFACE 

The  design  approach  is  determined  primarily  by  the  system  generating 
characteristics  and,  secondarily,  by  motoring  performance.  In  generating,  the 
system  must  convert  a raw  generator  output  which  varies  more  than  two  to  one  in 
voltage  and  frequency  to  regulated  28  Volts  DC.  As  a motoring  system,  the 
electronics  must  carry  out  both  commutation  and  motor  current  regulation. 

Several  possible  schemes  exist  to  accomplish  these  functions.  Three  trade- 
off areas  heavily  impact  the  system  design.  The  first  of  these  is  method  of 
commutation.  Define: 

Discrete  Commutation:  Motor  (or  generator)  commutation  by  sequentially 
connected  each  motor  phase  to  a fixed  voltage  source  such  that  the 
desired  direction  of  rotation  is  attained.  Often  called  "Square 
Wave  Drive". 

Continuous  Commutation:  Effecting  sinusoidal  current  drive  by  using 
pulse  width  modulation  (PWM)  techniques. 

Although  simple,  there  is  little  place  for  two  phase  discretely  commutated 
motors  in  high  efficiency  systems  due  to  the  third  harmonic  loss.  A three 
phase  approach,  however,  results  in  elimination  of  the  third  harmonic. 

A second  area  of  decision  is  Integral  versus  External  voltage  regulation; 

Integral  Regulation:  Use  of  pulse  width  modulation  to  combine  motor/ 
generator  commutation  with  current  regulation. 

External  Regulation:  Commutator  and  voltage  regulator  are  separate  entities. 

Integral  regulation  has  the  advantage  of  high  efficiency  since  most  driver 
losses  are  associated  with  the  voltage  drop  across  "ON"  semiconductors,  and 
this  approach  minimizes  the  number  of  in  series  semiconductor  switches  in  a 
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Fig.  A-l.  Performance  and  outline  drawing  of  brushless  dc  motor  generator. 


system.  It  should  be  emphasized  that  this  discussion  is  limited  to  digital 
drive  techniques  where  regulation  is  accomplished  by  pulse  width  modulation 
such  that  efficiency  is  much  higher  than  with  a linear  type  voltage  regulator. 
External  regulation  frees  the  commutator  to  be  optimized  for  the  commutation 
task,  such  as  with  very  high  motor  frequencies  where  the  PWM  frequency  gets  too 
high  and  switching  losses  become  significant. 

The  final  trade-off,  and  the  one  not  considered  in  the  proposal,  might  be 
called  High  Voltage  versus  Low  Voltage.  Refine: 

High  Voltage  System:  Motor /Generator  peak  back  emf  always  higher  than  the 
bus  voltage  (+28  VDC) . 

Low  Voltage  System:  Motor/Generator  back  emf  always  lower  than  the  bus 
voltage. 

SFS  has  done  considerable  work  with  the  latter  system.  As  will  be  shown, 
the  high  voltage  approach  offers  excellent  efficiencies.  A system  in  which  the 
peak  back  emf  may  vary  between  greater  than  and  smaller  than  the  bus  voltage 
presents  formidable  implementation  problems  in  the  electronic  driver. 

There  are  eight  combinations  of  these  three  parameters  that  a system  might 
have.  Most  will  be  rejected  by  the  following  arguments.  First,  consider 
Continuous  Commutation  with  the  other  two  choices.  Continuous  Commutation  with 
External  Regulation  is  redundant  since  the  commutation  process  itself  must 
regulate  the  bus  voltage  to  provide  sinusoidal  motor  current.  Continuous  Commu- 
tation in  a High  Voltage  System  is  possible,  but  requires  a very  complex  driver 
to  meet  the  combinations  of  motoring  or  generating,  instantaneous  bus  voltage 
higher  or  lower  than  instantaneous  back  emf,  and  instantaneous  motor  current 
positive  or  negative.  This  driver  is  not  only  complex,  additional  switch  drops 
are  introduced  which  lowers  its  efficiency.  Even  with  considerable  effort,  a 
driver  circuit  could  not  be  conceptualized  which  did  not  have  bidirectional 
power  switches  (switches  which  could,  under  control  of  a low  level  signal, 
either  block  voltages  of  both  polarities  or  pass  current  in  both  directions). 
Such  switches  can  be  built,  and  consist  of  two  transistors  and  two  diodes  each, 
with  up  to  2.9  volts  drop  across  the  switch.  Controlling  the  switches  is  also 
complex. 


94 


Because  of  the  additional  switch  drops,  a Low  Voltage,  Externally  Regulated 
system  cannot  compete  with  a Low  Voltage  system  with  Integral  Regulation.  A 
Discrete,  Integral,  High  Voltage  system  has  many  of  the  problems  of  a Continuous, 
Integral,  High  Voltage  System,  and  can  be  rejected.  A Discrete,  Integral,  Low 
Voltage  system  would  be  a viable  candidate;  a three  phase  version  using  fewer 
switches  than  its  two  phase,  Continuously  Commutated,  Integrally  Regulated,  Low 
Voltage  counterpart.  The  three  phase  design  would  be  slightly  less  efficient, 
however,  due  to  the  nonsinusoidal  current  waveform.  In  view  of  Sperry's  con- 
siderable experience  with  the  two  phase  approach,  this  type  of  discretely  com- 
mutated system  is  therefore  rejected. 

The  results  of  this  narrative  is  summarized  in  Table  A-2.  It  can  be  seen 
that  two  candidates  are  left: 

Type  4.  Discretely  Commutated,  Externally  Regulated,  High  Voltage  System 

Type  5.  ~:ly  Commutated,  Integrally  Regulated,  Low  Voltage  /stem 

The  latter's  opt.iiiu  .orm  factor  is  two  phase  with  "H"  bridges  for  each  phase. 
Sperry  experience  has  been  concentrated  in  that  area.  The  former  is  new  to 
Sperry. 

Figure  A-2  gives  a conceptual  schematic  for  the  new  approach.  Current 
sensing  techniques  are  not  ahown  and  would  have  to  be  worked  out.  Transistor/ 
rectifier  parts  Q1/CR1  through  A6/CR6  form  the  commutator  for  a three  phase 
delta  connected  motor.  As  stated  previously,  the  three  phase  conf iguration 
reduces  the  harmonic  content  of  the  discretely  commutated  drive  waveform.  The 
delta  format  is  chosen  to  minimize  commutator  requirements  in  that  the  delta 
does  not  require  one  of  the  transistors  of  a totem  pole  pair  to  turn  ON  just 
as  the  other  is  turning  JIT.  Many  "shunt  spike"  problems  are  therefore  <oided. 
The  number  of  motor  turns  is  also  maximized  with  the  delta  connection.  The 
six  diodes  full  wave  rectify  the  output  in  the  generation  mode.  The  motor 
interfaces  with  a high  voltage  Unregulated  Bus  which  is  roughly  equal  to  the 
peak  back  emf  of  the  motor. 

A pulse  width  modulated  switching  regulator  with  three  modes  of  operation 
is  formed  by  the  remaining  semiconductors.  A single  energy  storage  inductor. 
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LI,  is  used  in  the  regulation  operation.  When  in  the  generating  mode,  the 
system  Unregulated  Bus  is  stepped  down  by  the  PWM  action  of  07,  with  CRB  acting 
as  the  catch  diode.  When  motoring  in  the  normal  state  of  the  Unregulated  Bus 
being  higher  than  +28  volts,  a pulse  width  modulator  is  formed  by  08  and  CR7. 
Finally,  Q9  and  CR10  form  a third  pulse  width  modulator  v/hich  is  only  used  at 
system  start  up  when  the  motor  back  emf  is  less  than  the  bus  voltage.  Since  the 
Q9  modulator  is  only  required  at  start  up,  K1  can  be  used  to  increase  efficiency 
by  bypassing  that  regulator.  CR9  protects  09  in  case  K1  is  left  open  and  the 
system  returns  current  to  the  bus. 

The  Appendices  go  through  efficiency  calculations  for  each  of  the  two  can- 
didate systems.  A peak  Unregulated  Bus  voltage  of  100  volts  was  somewhat 
arbitrarily  selected  for  the  High  Voltage  system,  based  on  existing  semiconduc- 
tor technology  and  the  potential  high  reliability  application.  The  Low  Voltage 
approach  is  based  on  a back  emf  compatible  with  a 27  volt  minimum  bus. 

Table  A-3  presents  the  results  of  these  calculations.  The  range  of  efficien- 
cies for  each  mode  and  type  is  the  result  of  speed  dependencies.  A 2:1  speed 
range  C 7 5%  of  available  energy  extracted)  was  selected  for  this  comparison. 

The  very  attractive  generating  efficiency  of  the  High  Voltage  system  makes  it 
the  prime  candidate  for  implementation.  This  choice  is  considerably  strengthened 
by  the  very  high  pulse  width  modulation  frequency  required  with  a 50,000  RPM 
integrally  regulated  motor  and  by  an  initial  motor  design  which  revealed  dif- 
ficulty in  keeping  the  motor  back  emf  low  enough  for  the  Low  Voltage  approach. 

Although  not  shown  in  Fig.  A-2,  it  is  likely  that  additional  inductances 
will  have  to  be  added  to  the  commutation  circuit  to  control  semiconductor  cur- 
rents as  the  proposed  motor  has  negligible  winding  inductance.  Alternatively, 
it  might  be  possible  to  eliminate  Cl.  Since  the  winding  resistance  is  also 
virtually  negligible,  careful  circuit  design  will  be  necessary  to  prevent 
destructive  currents  from  flowing. 

The  motor /generator  selected  is  a three  phase,  delta  connected  device. 

The  motor /generator  will  operate  over  a 2.5:]  speed  range  (20K  or  50K  RPM). 

At  20,000  RPM,  the  device  as  a generator  must  have  sufficient  average  back  emf 
to  overcome  the  bus  voltage  plus  rectifier  drops.  Using  full  wave  rectification. 
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TABLE  A- 3 


DRIVER  EFFICIENCY  COMPARISON 


that  portion  of  a sine  wave  encompassing  60  degrees  each  side  of  peak  is  the 
active  waveform  to  be  rectified.  The  average  value  of  such  a waveform  is  0.8270 
times  the  peak  value.  Assuming  a 28  volt  bus  and  3 volts  of  switch/wiring  drop, 
the  back  emf  at  minimum  operating  speed  must  then  be  at  least: 

,,  28+3  ot  r , . 

V . = ■ , ,n-,n  = 37.5  volts  peak 

min  .8270  v 

As  an  alternative  to  an  auxiliary  start-up  circuit  to  accommodate  the  condi- 
tion of  motor  peak  back  emf  less  than  the  bus  voltage,  start-up  could  be  accom- 
plished by  wave  width  modulating  the  commutator  to  limit  the  current.  Efficiency 
in  this  mode  is  unimportant  since  it  is  only  done  once  during  the  flight.  This 
could  complicate  the  low  level  circuitry  somewhat,  but  eliminates  Q9,  GR9,  CR10, 
and  K1  from  Fig.  A- 2. 

Some  time  was  spent  looking  at  wave  width  modulation  techniques  for  the 
three  phase  concept.  By  using  sinusoidal  sense  winding  outputs,  it  was  con- 
cluded that  speed  independent,  wave  width  modulation  could  be  accomplished 
which  was  symmetrical,  eliminating  the  third  harmonic  under  all  conditions. 

A. 5 MOTOR/GENERATOR  DESIGN 

This  project  presents  a challenge  in  two  distinct  areas — electromagnetic 
and  mechanical.  In  arriving  at  a final  design,  both  areas  must  be  considered. 

The  next  section  will  present  a development  of  the  basic  equations  used  in  the 
design  of  this  motor/generator  followed  by  a section  on  the  final  design. 

A . 5 . A . Development  of  Models 

The  decision  to  use  an  external  inductor  as  part  of  the  commu- 
tation/regulation electronics  makes  modeling  the  inductance  of  less  importance 
than  anticipated  in  the  proposal.  The  large  air  gap,  absence  of  iron  in  the 
armature,  and  low  permeability  of  the  Samarium  cobalt  permanent  material 
inherently  result  in  very  low  inductance  which  can  be  neglected  in  the  motor/ 
generator  analysis.  Thus  the  primary  parameters  needed  to  establish  the  per- 
formance are  the  winding  resistance  and  the  generated  back  emf.  The  resistance 
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calculation  is  simply  a matter  of  determining  the  length  of  the  conductor  from 
the  geometry  of  the  armature  and  will  not  be  detailed.  The  back  emf  was  cal- 
culated as  follows: 

The  basic  equation  for  voltage  due  to  a conductor  moving  in  a magnetic 
field  is: 

e = B IV  x 10~8  (1) 

where  e = The  generated  voltage  in  volts. 

2 

B = Magnetic  flux  density  in  lines/in  . 

1 = Effective  length  of  conductor  in  B field  in  inches. 

V = Velocity  of  conductor  relative  to  field  in/sec. 

the  velocity  due  to  rotation  is: 

V = wr  (2) 

w = Angular  velocity  in  rad/sec. 

r = Radius  in  inches. 

substituting  (2)  into  (1)  and  rearranging  - 

- = B 1 r x 10-8  (3) 

w 

for  multiple  conductors 

Kj-^Blr  N/eff  x 1(T8  (A) 

For  our  case,  we  need  only  to  determine  the  flux  density  in  the  air  gap  and 
the  effective  turns.  All  basic  discussions  of  permanent  magnetic  circuits  show 
that  minimum  magnetic  material  1b  required  to  obtain  a specific  air  gap  flux 
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density  if  the  magnet  is  operated  with  a shear  slope  to  give  maximum  energy 
product  (B  x H) . We  have  found  at  Sperry  that  the  minimum  total  weight  for 
! conventional  motors  and  generator  is  usually  obtained  at  somewhat  higher  shear 

slopes  and  never  at  lower.  For  this  design,  using  an  ironless  armature  result- 
ing in  a very  large  magnetic  gap,  and  requiring  substantial  sections  of  non- 
magnetic material  to  provide  mechanical  strength,  it  is  not  practical  to  utilize 
sufficient  length  of  permanent  magnet  material  to  obtain  a shear  slope  signifi- 
cantly greater  than  1.0  which  is  the  maximum  energy  point  for  Samarium  cobalt 
and  this  was  selected  for  this  design. 

The  magnet  flux  density  can  be  determined  from  the  manufacturer  supplied 
demagnetization  curve  for  a B/H  of  1 to  be  4K  Gauss.  Some  suppliers  provide  a 
curve  for  a "reluctance  factor"  to  reduce  the  flux  and  obtain  the  useful  flux. 
This  method  or  a leakage  factor,  usually  imperically  or  intuitively  obtained, 
must  be  used  to  obtain  the  "useful"  or  air  gap  flux  density.  Either  method 
results  in  a factor  of  0.8  or  an  air  gap  flux  density  of  4.0  x 0.8  = 3.2  Kilo 
Gauss  or  20  Kilolines  per  square  inch  magnetic  flux  density  in  the  gap. 

The  effective  number  of  turns  is  obtained  from  the  final  magnetic  con- 
figuration by  determining  the  maximum  number  of  turns  which  will  be  located  in 
the  B-field.  Using  this  number  of  turns  will  give  the  peak  back  emf. 

The  major  mechanical  design  problem  was  to  contain  the  rotating  permanent 
magnets  at  speeds  of  50,000  RPM.  The  analysis  of  the  stresses  in  the  rotor  were 
conducted  using  standard  methods,  primarily  from  Roark's  book,  "Formulas  for 
Stress  and  Strain". 

A.5.B.  Final  Design  Parameters 

1.  Electromagnetic  Design  and  Performance  Calculations 

The  electronic  interface  analysis  (Section  A. 4)  defined  a 
three-phase  motor/generator  with  a voltage  of  37.5  volts  peak  minimum  at 
20,000  RPM  or  a back  emf  constant  of  - 

h “ 20,00F~x".1047  “ -°18  Volts/Rad/Sec  minimum 
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For  maximum  efficiency,  the  winding  resistance  should  be  minimum. 
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Analysis  of  Generator  Losses 

The  major  losses  in  a conventional  permanent  magnet  generator  are: 

2 

Stator  or  Armature  I R losses 

Fundamental  Frequency  Iron  losses 

(Eddy  Current  and  Hysteresis) 

High  Frequency  Iron  losses 

For  the  ironless  armature  design,  the  fundamental  iron  loss  is  nonexistent.  The 
high  frequency  iron  loss  is  caused  by  pulsations  of  the  flux  as  the  magnet  rotates 
past  the  slot  openings.  This  is  also  nonexistent.  At  the  high  frequency  encoun- 
tered at  high  speed,  consideration  must  be  given  to  increased  resistance  of  the 
conductors  due  to  skin  effect.  Also  since  the  conductors  are  located  directly 

in  the  flux  field,  the  effect  of  eddy  current  losses  must  be  evaluated. 

2 

The  losses  to  be  considered  then  are  the  armature  I R losses  and  eddy  cur- 
rent losses  in  the  conductors.  The  original  concept  for  the  winding  was  a delta 
connected  winding  using  five  active  poles  for  each  redundant  half  and  one  pole 
for  each  sense  winding.  After  some  preliminary  measurements  of  air  gap  flux 
density  showed  approximately  10 7,  variation  from  the  magnet  surface  to  the  return 
path,  it  was  concluded  that  the  windings  should  be  connected  in  wye  to  elimi- 
nate the  possibility  of  losses  due  to  circulating  currents. 

Attempts  to  calculate  the  eddy  current  losses  due  to  the  changing  flux 
density  proved  fruitless.  By  using  several  "reasonable"  assumptions,  the 
calculated  losses  varied  from  neglectible  to  greater  than  rated  output.  A 
test  was  conducted  to  determine  if  these  losses  were  significant.  A toroidal 
core  of  laminated  nickel  iron  was  fabricated  as  follows: 


Outside  Diameter 

Inside  Diameter 

Length 

Slot 

Turns 

Resistance 

Inductance 


6 inches 
4-7/8  inches 
11/16  inches 
1/16  inch 
100 

0.4  ohm 
3.92  mH 
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Tests  were  conducted  to  determine  the  flux  density  in  the  Rap  versus  current  as 
follows : 


Input  Current  Gap  Flux  Density 

0.25  270  Gauss 

0.50  Amp  553  Gauss 

0.75  Amp  818  Gauss 

Tests  were  then  made  to  measure  the  input  power  when  excited  at  2000  Hertz  and 

5000  Hertz  with  nothing  in  the  gap,  with  one  of  the  coils  wound  with  Litz  wire, 

with  a coil  wound  with  solid  wire,  and  with  a sheet  of  aluminum  foil.  The  cur- 
rent was  set  at  .65  amp  rms  to  give  a gap  density  of  1KC-  peak.  The  increase 

in  power  as  the  various  samples  were  placed  in  the  air  gap  is  tabulated  below: 


IN  GAP 

INCREASED  POWER  AT  2K  HZ 

INCREASED  POWER  AT  5K  HZ 

Nothing 

0 

0 

Coil  of  Litz  Wire 

<.03 

<.03 

Coil  of  Solid  Wire 

<.03 

<.03 

.002  Th.  Aluminum  Foil 

>2.0 

>4.0 

.005  Th.  Aluminum  Foil 

>4.0 

>8.0 

In  the  above  test,  no  increase  in  power  could  be  measured  with  the  generator 
coils,  either  using  solid  wire  or  Litz  wire.  The  resolution  of  the  meter  was 
approximately  .03  watts.  As  a worst  case,  then,  the  eddy  current  losses  will 
be: 


Max.  Losses 


( _B_  ) 

( 1.0  ) 


2 

x Number  of  Coils  x .03 


( a 2 r 

Max.  Losses  ■ ^ y— ^ x 12  x .03  * 3.7  Watt 


The  increase  in  resistance  due  to  skin  effect  for  solid  wire  is: 
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SPEED 


AC  RESISTANCE 


20.000  1.2  x Rdc 

50.000  2.25  x Rdc 

With  forty  strand  I.Itz  wire,  the  increase  in  resistance  is  neglectible. 

A summary  of  the  calculated  parameters  for  the  generator  is  tabulated  in 
Table  A-4  using  tested  values  for  selected  parameters.  This  table  shows  that 
the  generator  efficiency  is  92.9%  at  20,000  and  97.6%  at  50,000  RPM  using  the 
worst  case  eddy  current  losses.  These  calculations  are  based  on  one  of  the 
two  redundance  windings. 

Performance  as  a motor  is  summarized  in  Table  A-5  using  the  same  basic 
parameters.  This  table  shows  that  the  motor  efficiency  is  comparable  to  the 
generator  under  the  condition  of  250  watts  constant  output  power.  Tf  it  is 
assumed  that  the  duty  cycle  specified  in  "Specification  for  Flywheel  Motor- 
Generator  Design"  is  used  and  the  average  torque  for  constant  acceleration  is 
used,  the  efficiency  drops  to  84%  at  20,000  RPM,  93%  at  50,000  RPM,  or  88.6% 
average  efficiency.  The  primary  loss  in  this  case  is  the  eddy  current  loss 
of  3.7  watts  which  is  a worst  case  number  and  is  likely  high. 

2.  Mechanical  Design  and  Stress  Calculations 

The  motor  generator  consists  of  three  separate  pieces — the 
armature  or  stator,  the  permanent  magnet  rotor,  and  the  return  path.  See 
Outline  Drawing  in  Fig.  A-l.  The  construction  of  each  of  these  will  be  dis- 
cussed separately  below. 

a.  The  Armature 

The  armature  is  fabricated  using  three  coil  forms 
made  of  epoxy  impregnated  fiberglass  to  hold  and  provide  proper  spacing  of 
the  twelve  wound  coils.  The  coils  are  wound  on  a special  form  to  give  the 
same  shape  as  the  coil  form  and  are  impregnated  prior  to  removing  from  the 
winding  form.  After  the  coils  are  assembled  and  bonded  to  the  coil  forms,  the 
three  stages  are  located  in  an  aluminum  bracket  and  the  coils  are  connected 
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TABLE  A-4 


GENERATOR  PERFORMANCE 


PARAMETER 

Winding  Resistance 
Back  eraf  - Kg 
Inductance 


UNITS 


Volts/Rad/Sec 

Microhenries 


CALCULATED  VALUE 


At  20,000  RPM 


Generator  Voltage 
Average  Voltage 

Current  at  250  Watts 

2 

I R Losses 
Eddy  Current  Losses 
Total  Losses 
Efficiency 

At  50,000  RPM 

Generated  Voltage 

Average  Voltage 

Current  at  250  Watts 
2 

I R.  Losses 
Eddy  Current  Losses 
Total  Losses 
Efficiency 


Volts  Peak 
Volts  Peak 
Amps  DC 
Watts 
Watts 
Watts 


Volts  Peak 
Volts  DC 


Watts 

Watts 

Watts 


120.4 


Average  Efficiency 


* ( 


TABLE  A- 5 


MOTOR 

PERFORMANCE 

PARAMETER 

UNITS 

CALCULATED  VALUE 

Winding  Resistance 

Ohms 

0.39 

Back  emf  - 

Volts/Rad/Sec 

.023 

Inductance 

Micro  Ohms 

15 

Torque  Constant 

Oz.  In/Amp 

3.26 

Average  Torque  for  10  Hr.  Run-up 

Oz.  In. 

1.35 

At  20,000  RPM,  1.35  Oz.-In. 

Generated  Voltage 

Volts  Peak 

48.1 

Input  Current 

Amps 

.414 

2 

I R Losses 

Watts 

.067 

Eddy  Current  Losses 

Watts 

3.7 

Total  Losses 

Watts 

3.77 

Output  Power 

Watts 

20 

Efficiency 

% 

84.1 

At  20,000  RPM,  250  Watts 

Generated  Voltage 

Volts  Peak 

48.1 

Input  Current 

Amps 

6.28 

2 

I R Losses 

Watts 

15.4 

Eddy  Current  Losses 

Watts 

3.7 

Total  Losses 

Watts 

19.1 

Output  Torque 

Oz.  In. 

16.9 

Efficiency 

% 

92.9 

At  50,000  RPM,  1.35  Oz.-In. 

Generated  Voltage 

Volts  Peak 

120.4 

Input  Current 

Amps 

.414 

I^R  Losses 

Watts 

.067 

Eddy  Current  Losses 

Watts 

3.7 

Total  Losses 

Watts 

3.77 
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TABLE  A- 5 (Continued) 


PARAMETER 

UNITS 

CALCULATED  VALUE 

Output  Power 

Watts 

50 

Efficiency 

®/ 

/© 

93 

At  50,000  RPM,  250  Watts 

Generated  Voltage 

Volts  O-Peak 

120.4 

Input  Current 

Amps 

2.52 

2 

I R Losses 

Watts 

2.48 

Eddy  Current  Losses 

Watts 

3.7 

Total  Losses 

Watts 

6.18 

Output  Torque 

Oz.-In. 

6.77 

Efficiency 

% 

97.6 

Average  Efficiency  at  1.35  Oz-In. 

’A 

88.6 

Average  Efficiency  at  250  Watts  Out 

Z 

94.5 

Sense  Coil 

Resistance 

Ohms 

.05 

Back  emf  Constant  - 

Volts/Rad/Sec 

.005 

using  press  fit  terminals  around  the  outer  land  of  the  bracket.  The  jumpers 
and  connectors  are  insulated  and  protected  with  a polyurethane  coating.  The 
coil  assembly  is  then  potted  using  a heavily  filled  epoxy. 

b . The  Return  Path 

The  return  path  is  machined  from  a high  strength  stain- 
less steel  and  heat  treated  for  maximum  strength. 

c . The  Permanent  Magnet  Rotor 

The  permanent  magnet  rotor  presented  the  primary  chal- 
lenge to  the  mechanical  design.  After  several  false  starts,  the  final  design 
was  selected  to  consist  of  a high  strength  titanium  ribbed  structure  to  contain 
the  magnetics  on  five  sides,  a magnetic  plate  containing  the  magnets  on  the 
final  side  and  providing  a return  path,  and  a final  titanium  ring  with  a shrink 
fit  to  share  the  stress  load.  The  sketch  in  Fig.  A-3  shows  the  general  con- 
struction of  this  member.  Some  of  the  key  design  points  are: 

Outer  surface  of  magnet  flat  to  distribute  load 

Magnets  bonded  at  OD  only  (i.e..  All  stresses  are  compressive) 

Return  path  supported  at  ID  only  (i.e.,  does  not  induce  any 
stress  load  in  Ti  housing) . 

Magnet  is  contained  on  all  sides. 

Appendix  3 contains  the  Rotor  Stress  Analysis.  A summary  of  the  stresses 
are  shown  in  Table  A- 6. 


TABLE  A-6 
STRESS  SUMMARY 

MAXIMUM  STRESS  (KPS1) 

17.7  (Compresses) 
94.2  (Tension) 

50  (Tension) 

42  (Tension) 


SECTION 

Magnetic 
Outer  Ring 
Ribs 

Inner  Hub 
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WEIGHT 


The  calculated  weight  for  the  motor  generator  Is  four  pounds  with  a break- 
down as  follows: 


WEIGHT 

1.0 

.42 

.42 

.87 

2.71 


Return  Path 

.71 

Wire 

.20 

Coil  Form 

.15 

Coll  Bracket 

.16 

Potting  and  Miscellaneous 

.1 

Total  Stator 

.61 

TOTAL 

4.03 

ITEM 

Rotor  Shell 
Magnet  Plate 
Magnets 
Outer  Ring 

Total  Magnet  Rotor 
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A. 6 COMMENTS  ON  SAMPLE  FABRICATION 


In  general,  the  sample  worked  out  well  as  designed  with  only  minor  cor- 
rections being  required.  The  major  item  which  will  require  changing,  if 
additional  units  are  built,  is  to  increase  the  diameter  of  the  armature  flange 
to  allow  more  room  for  the  wire  terminations.  Difficulty  was  encountered  in 
trying  to  contain  this  wire  and  it  was  necessary  to  omit  sleeving  from  these 
jumpers  on  the  breadboard  unit. 

Material  substitutions  were  made  on  this  unit  to  use  material  in  stock  at 
Sperry.  Thus,  the  rotor  housing  is  fabricated  from  30A  Stainless  Steel  rather 
than  Titanium  and  the  magnetic  plates  are  A16  steel  rather  than  AAOC.  These 
substitutions  would  prevent  any  high  speed  test. 

The  test  performed  on  the  sample  at  Sperry  consisted  of  measuring: 

Resistance 
Inductance 
Back  emf  constant 
Wave  shape 

A summary  of  these  data  is  presented  in  Table  A-7  and  Figs.  A-A  and  A-5. 

Photos  of  the  engineering  unit  are  shown  in  Figs.  A-5  through  A-8. 


TABLE  A- 7 


TEST  DATA  ON  BREADBOARD  UNIT 


TERMINALS 

RESISTANCE 

(Ohms) 

INDUCTANCE 

(mH) 

•s 

(V/Rad/Sec) 

Blk  (18)  - Red  (5) 

.396 

.0344 

.0239 

Red  (5)  - Blu  (4) 

.397 

.0358 

.0243 

Blu  (4)  - Blk  (18) 

.399 

.0347 

.0244 

Blk  (36)  - Red  (23) 

.388 

.0353 

.0247 

Blu  (22)  - Blk  (36) 

.385 

.0338 

.0246 

Wht  (2)  - Yel  (3) 

.0853 

. 0094 

.00499 

Yel  (3)  - Grn  (1) 

.0841 

.0089 

.00490 

Grn  (1)  - Wht  (2) 

.0858 

.009 

.00488 

Wht  (20)  - Yel  (21) 

.090 

.009 

.00499 

Yel  (21)  - Grn  (19) 

.085 

.009 

.00495 

Grn  (19)  - Wht  (20) 

.087 

.0091 

.00488 
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5 V / cm 

OUTPUT  = 17  V P-P 


1 V / cm 

VARIATION  < 0.2  V 
<2.3  PERCENT 


Fig.  A-5.  Output  uniformity. 


r 
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A . 7 Two-Phase  Cont i nuously  Commutated  Efficiency  Calculations 
TWO  PHASE  BACK  EMF  < 28V  EFFICIENCY 
A.  Motoring 


Fig.  A-9.  Motoring  diagram. 


Sinusoidal  current  is  pumped  through  the  motor  by  holding  S4(S2)  on  and  Pul  Be 

Width  Modulating  SI (S3)  (Parentheses  mode  is  negative  current  flow;  S1-S4  is 

defined  as  positive  current  flow.)  Consider  the  positive  half  cycle.  Assume 

that  the  average  of  switch  drop  and  catch  diode  drop  is  1.5V.  Therefore,  the 

average  switching  voltage  drop  in  always  twice  that,  or  3.0  Volts. 

The  power  delivered  to  the  motor  which  is  converted  to  mechanical  power 

is  i e.  , where  e,  is  the  motor  back  emf: 
mb  b 


eb  - KLj^O)  SIN  wT 

Where  - Back  emf  constant  in  volts/rad/sec. 

oj  - Motor  rotational  velocity  in  rad/sec. 

2 

The  motor  copper  loss  is  i R . 

mm 

And  the  switch  voltage  drop  loss  is  1 V^ 
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Where  = Average  switch  drop 

*=  3.0  Volts  as  discussed  above. 

It  is  further  assumed  that  the  driver  loss  due  to  switching  transients  is 
negligible.  This  is  a warranted  assumption.  Without  going  through  the  deriva- 
tion, the  loss  in  any  one  switch  is  approximately 


Where  V = Bus  Voltage 

n 


i - Average  Motor  Current 
m 


t *=  Switch  rise  (or  fall)  time 

T “ PWM  period 


It  can  easily  be  shown  that  i ■ — I , 

m tt  pk 

Where  I . is  the  peak  motor  current, 
pk 


i - I . SIN  (joT 
m pk 

tt  1 V 

- COS  WT 
TT 
m 


T * 

■-?  I 


SIN  cjT  dwt 


- - 1 
it  pk 


Note  that  i is  selected  to  be  in  phase  with  the  back  emf.  This  results  in 
in 

maximum  efficiency. 

Taking  both  SIN  and  COS  phases  into  account: 

6 1.  g i 

Pout  ■ Mechanical  Output  Power  » b m + b m 

s s c c 


Where  eb 

s 

i 

m 

s 


Back  EMF  for  SIN  winding 
Motor  Current  in  SIN  winding 


eb^  - Back  EMF  for  COS  winding 
*mc  ■ Motor  Current  in  COS  winding 
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Pout  = (Kh  OJ  SIN  uT)  (Ipk  SIN  0jT)  + (1^  w COS  wT)  (IpR  COS  U)T 


Kh  uj  Ipk  (SIN2  ujT  + COS2  tuT)  = (j  I 

Motor  Losses  (ignoring  iron  loss)  = Pml 


Pm  I = i R 
m m 


Pml 


R " 

~5  / 


TT 


i 2 dwt  - — I . 2 / SIN2  wT  dwt 

m TT  pk  J 


R I , 
tn  pk 

IT 


[1/2  o>T  - 1/4  SIN  2 coT] 


R 1 . R I . 

m pk  . ii  . m pk 

tt  1 2 J 2 


(Per  Phase) 


And  driver  loss  = P,_  is: 

dl 


dl 


dl 


^ VD 


/ 1 

u 


2 V 1 

I k SIN  toT  dwt  = ^ ^ (Per  Phase) 


Bus  Power 


P = Pout  + Pml  + P,, 

dl 


V 1 - + R 1 ^ + ^ v-  ± 1 


2 + i 

pk  n ‘pk  tt  D 

I , [K,  cj  + R I , + 1 VJ 
pk  d m pk  tt  D 


pk’ 


And  Driver  Efficiency 


nd  = x ioo% 


N, 


KBto  + R I , 

B m pk 

D K_u>  + R I . + - V„ 
B m pk  tt  D 
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For  further  simplification,  assume  Pml  » 0 


Then 


S' 


to 


V + ¥ VD 


With  the  assumption  that  switch  and  diode  drags  are  equal  (in  reality, 
diode  drop  is  ~ IV  and  switch  drop  is  1.9V),  the  same  equation  applies  to 
generation.  Referring  back  to  Fig.  1,  generation  is  accomplished  by  maintain- 
ing — 


i « I , SIN  toT 
m pk 

This  is  done  by  opening  all  switches  except  S4(S2),  which  is  pulse  width 
modulated.  On  closing,  the  motor  is  effectively  shorted,  and  im  builds  up.  On 
opening,  the  energy  stored  in  provides  an  output  to  the  bus  via  CR3  (CR1), 
even  though  the  back  emf  is  less  than  the  bus  voltage.  Generating,  in  general, 
will  be  slightly  more  efficient  than  motoring  because  of  the  substitution  of  a 
diode  drop  for  a switch  (transistor)  drop  in  the  former  case.  ( 2%  higher) 

So,  let's  design  a motor.  The  back  emf  must  always  be  less  than  the  bus 
voltage  minus  switch  drops  so  that  motoring  may  be  accomplished. 


Set 

(min) 

27  Volts 

VD 

- 

3 Volts 

VF 

- 

Wiring,  Filter  Drop  • IV  N.C. 

Then 

eb(max) 

■ 

KgtoW  max.  = 27-3-1  “ 23  Volts 

to 

If  to 

min 

- 

max 

2 

Then 

min 

- 

11.5  Volts 

And: 
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KRu  min 


K <jj  max 
B 


x 100% 


V min  + ¥ vo 


, x 100%  < N < 
4 „ — I)  — 


V n^x  + - VD 


— x 100%  < N < x 100% 

4 „ — D — , A 


11.5  + - .3 

TT 


23  + - .3 

TT 


75%  1 ND  1 86% 


Readdressing  switching  time  losses: 

Nc  = Switching  time  efficiency 

— x 100% 


kb“ 


2V 


S' 


w + 


B .2  tr 


TT.  6 


If  T 

= 50  M sec 

tr 

= 1 M sec 

Then, 

at  min;  N<,  = 

11.5 

n-5  + ff  • 

2 

50 

X 100% 


99% 
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A. 8 Three-Phase,  Discretely  Commutated  Efficiency  Calculations 


THREE  PHASE  EFFICIENCY 
Def Inltlons 

V = Bus  voltage,  nominalcy  +28V 

D 

= Unregulated,  or  Motor/Generator  bus  voltage. 

ICg  = Back  EMF  constant,  V/Rad/Sec. 

W = Rotational  Velocity,  Rads/Sec. 

- Voltage  drop  across  transistor  switch,  < 1.9V. 

Vp  = Voltage  drop  across  diode,  < 1.0V. 

I„  “ +28V  Bus  current. 

D 

I„  = Motor  (unregulated  bus)  current. 

N - Efficiency,  with  subscripts  delineating  commutator,  regulator, 

generating,  and  motoring  modes,  i.e.:  N indicating  efficiency 
of  the  commutation  in  the  generating  mode? 
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THREE  PHASE  EFFICIENCY  (Neglecting  external  inductor  losses) 
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Now, 


vu  “ ^ w 

Realizing  the  state-of-the-art  in  semiconductor  technology,  combined  with  the 

potential  high-reliability  application,  set  Vy  (Max  ) “ 100  volts.  If  ajmln  " 

1/2  oj  , then  50  < V,T  < 100V.  For  w , : 
max  — U — min 

<,  2 1 9, ri_ t1-’  «-f§> + 1-° <n>]n 

»n  ■ 1 2T ^Jj  *10<“ 


(.924)  (1  - 


28 

.836  + .56 


28 


) x 100% 


- (.924)  (.950)  x 100%  - 88% 


For  0) 


max 


„ 2.19.  , <1-9<1-TO)  + 1-0(llo) 


100 


28 


x 100% 


- (.962)  (1  ~ 1v36281'  ,28>  x 100% 


- (.962)  (.941)  x 100%  - 91% 


B . Generating 

1.  Commutation  Efficiency,  N 

» - 

Cg  \ vu  / 


eg 


2.  Regulator  Efficiency,  N 


r8 


Q7  and  CR8  form  the  pulse  width  modulator  in  this  mode. 

vB 

Q7  is  on  — of  the  time. 

U 

vB 

CR8  conducts  1 — tt-  of  the  time. 

U 
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Losses  are 


[VT^V0 


Total  output  power  i6  Vglg>  80 


rg 


V 


1 - 


V VB  + V I 1 — 

LT\  v» 


B 


x 100% 


and  combined  efficiency,  N„  is: 

o 


2V_ 


1 - 


1 - 


X 100% 


For  w 


min' 


u 50 


(.960)  (1  - •■-)  x 100% 


- (.960)  (.946)  x 100%  - 91% 


For  oj  ; 
max 


2.i  r [i>9  loo + 1,( 

100  23 


(0.98)  (1  - ,53-22g— 72)  x 100% 


(1 


0 <1  - T5o>l 


- (.98)  (.955)  x 100%  - 94% 


100% 


x 100% 
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A.  9 Rotor  Stress  Calculations 


STRESSES 


Rotational  Stress 


O.D. 


I.D. 


A 

R 

W 

S 


S 


5.3 


A. 875 

.175 

2.5AA 


5235  rad/sec. 


AR^2 
386. A 


(page  360) 


.175  x 2.5AA2  x 52352 
386. A 


80.3  K 


psi 


Mass  of  Magnet 


W 

F 


.85  + . A5 
2 

My  W2  - 


x .75  x 

Wy  w2 
g 


.25  x 


.0353  x 1.5  x 52352 
386. A 


29  - .0353  lb.  each 


3.755  Klb.  each 


Case  6,  Page  217 


S 

max 


3W 

? 

Awt 


J 


V 


17.67  K psl 


m = 

1/V  - 

3 

Totals 

St  = 

80.3  + 

9.5  = 89.8 

Sy  = 

28.46 

S = 

|^9.82 

+ 28. 462 

S = 

94.2 

W 


Unit  Load  = 


3.755  x 10 
.85  x .25 


t = .385 


r max 


3 x 17.67 
4tt  x .385^ 


28.46  K psi 


st  = 


3W 


r max  28.46 


47Tt 


In  Ribs 

F 3.755  K psl 

° ” A .3  x .25 

0 = 50.06 


Inner  Hub 

Rotational 

0. D.  - 3.0 

1. D.  - 2.0 

L - .625 
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386.11 


.175  x 1.252  x 52352 


386.4 


19.39  Kpsl 


Case  9,  Page  174 


Max.  Tension  - 1/2  W COT  0 


- 1/2  x 


3.755  x 
tan  15° 


7.01  Klb. 


7.0  Klb. 
.625  x .5 


- 22.4  Kpsi 


Total  22.4  H 19.39  - 41.8  O.K. 


STRESS  ON  MAGNET 


F - 3.755  K each 


A - .850  x .250  - .2125 


o ■ F/A  ■ 17.7  Kpsl.  compressure 

VS  42 


NOTE:  All  page  numbers  and  case  numbers  refer  to:  Formulas  for  Stress  and 
Strain,  R.  J.  Roark  (McGrav  Hill  Publishing  Company)  Fourth  Edition, 
1965. 
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A . 10  Addendum — Final  Report 
MAGNETIC  CIRCUIT 

The  magnetic  circuit  is  designed  to  utilize  the  high  coercivlty  of  Samarium 
cobalt.  As  a result,  manufacturing  tolerances  and  rotor-stator  clearances  do 
not  appreciably  affect  performance.  Figure  A-10  shows  the  basic  magnetic  circuit 
including  the  dimensions  affecting  the  magnet  operating  point. 


6 H9ie 


rotor 

SHELL 


qm  MECHANICAL 
AIR  GAP 


ROTOR 

PLATE 

\ 

X 


STATOR 

l 


MAGNE  TIC 
PLATE 


0.250  t 0.001 
0 015  ± 0 002 


qm 


0 250  i 0 001 


Fig.  A-10.  Magnetic  circuit. 
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The  effective  magnetic  air  gap  is  the  sum  of  the 
magnet,  the  stator  thickness,  and  2 gm.  For  purposes 
cal  air  gap  is  assumed  to  be  .020  ± .005.  From  worst 
netic  air  gap  is: 

Min.  Gap  ■ 0.292  in. 

Max.  Gap  ■ 0.318  in. 

Thus,  the  extremes  of  magnet  operation  are: 


Min.  B/H 
Max.  B/H 


Min.  Magnet  Length 
Max.  Gap  x Leakage  Factor 

Max.  Magnet  Length 
Min.  Gap  x Leakage  Factor 


rotor  shell  covering  the 
of  analysis,  the  mechani- 
case  tolerances,  the  mag- 


.249 

.318 

.251 

.292 


x . 80 


x . 80 


.98 


1.075 


From  the  normal  demagnetization  curve  for  Samarium  cobalt,  the  operating  flux 
density  is: 

B...  ■ 3.95  Kilogauss 

fiin  • 

B„  *4.15  Kilogauss 

Max. 

This  variation  is  +3.75%,  -1.25%  from  the  design  nominal  of  4.0  Kilogauss.  Con- 
siderable latitude  is  thus  permissible  in  establishing  the  mechanical  air  gap 
without  significantly  affecting  performance. 

Another  factor  which  contributes  to  a conservative  design  is  that  Samarium 
cobalt  magnets  manufactured  recently  are  exceeding  the  normal  demagnetization 
curve  used  (8  Kilogauss  to  8 Kilo-oersteds). 

WORST  CASE  EFFICIENCY 

Efficiency  of  the  unit  as  a generator  is  controlled  by  three  factors: 

1.  Generated  Back  EMF, 

2.  Winding  Resistance,  and 

3.  Rotational  Losses. 

Under  the  operating  conditions  of  20,000  to  50,000  RPM,  the  variation  in 
generated  back  EMF  is  controlled  by  the  air  gap  flux  density.  Variation  in 
density  due  to  manufacturing  tolerances  and  air  gap  dimensions  was  shown  to 
be  +3.75%,  -1.25%.  The  other  variable  is  magnet  fabrication  consistency. 
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This  is  typically  >52.  From  Sperry's  experience  with  Samarium  cobalt,  total 
unit-to-unit  variation  in  back  EMF  does  not  exceed  1 1 02.  This  factor  is  used 
for  worst  case  analysis  and  is  considered  quite  conservative. 

Winding  resistance  variation  is  primarily  due  to  wire  fabrication  toler- 
ances but  is  estimated  to  be  >10 7.  to  allow  for  winding  technique  and  jumper/ 
lead  wire  routing. 

Rotational  losses  (eddy  current  losses)  were  determined  as  a worst  case 
limit  in  the  final  report  and  are  projected  to  be  3.7  watts  maximum. 


PARAMETER 

NOMINAL 

WORST  CASF. 

Kg  = Back  FMF  (Volts/Kad/Sec. ) 

.023 

.0207 

R = Resistance  (Ohms) 

.39 

.429 

W = Eddy  Current  Losses  (Watts) 

ti 

— 

3.7 

At  20,000  RPM 

The  worst  case  voltage  at  20,000  RPM  (2094  Rad/Sec)  is: 
V = K„W  = .0207  (2094)  = 43.35  Volts  Peak 

D 

and  the  average  voltage  is  35.86  volts  peak. 

The  current  at  250  watts  is  then — 


250 


6.97  Amps. 


35.86 

Using  the  worst  case  resistance  of  .429  ohm,  the  I2R  loss  is  20.85  watts. 
Thus  the  worst  case  efficiency  at  20,000  RPM  is — 


Output  Power  x 100 
2 

Output  Power  + I R + W 

h. 


250  x 100 


250  + 20.85  + 3.7 


91.067 


At  50,000  RPM 

Using  the  same  calculation  method  as  for  20,000  RPM  results  in  the  follow- 
ing worst  case  parameters. 

Average  Voltage  = 89.65  Volts 

Current  at  250  Watts  **  2.789  Amps 
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2 

I R Losses  ■ 3.34  Watts 

Worst  Case  Efficiency  » 97.26 % 

The  average  efficiency  over  the  operating  range  is  then  95. 2Z  worst  case. 


134 


ACKNOWLEDGEMENTS 


The  author  is  grateful  to  Cornelius  M.  Sullivan  for  the  inspired  mechan- 
ical engineering  effort  resulting  in  the  magnetic  hearing  realization.  Also, 
thanks  are  due  to  Armand  Bouvier  whose  persistent  design  work  kept  up  with 
ever-changing  concepts;  and  to  the  support  of  the  Div.  6 Machine  Shop  under 
the  direction  of  Frank  Weibel  and  William  Fielding,  who  produced  excellent 
hardware  in  a timely  manner.  Thanks  also  to  Steve  Britten  for  the  magnetic 
finite  element  program  development,  and  to  David  Hyland  for  analysis  of  the 
anisotropic  rotors. 


UNCLASSIFIP.U 

’ t ■ jyRI  W rL^tHTTh,*  TION  OF  THIS  PAGE  (then  Onto  ( nlrrrd) 

(l  7 ) ( \C1  J REPORT  DOCUMENTATION  PAGE 


TITLE  (anil  Subtitle) 


2.  GOVT  ACCESSION  NO 


Flywheel  Components  for  Satellite  Applications  • 


KE.AI)  INSTKIJf  TIONS 

BKEOkt.  f OMPI.KTINO  FORM 

3.  RECIPIENT'S  CATALOG  NUMBER 


t ip  E OF  REPORT  & PERIOfi.  COVERED 


/Cf  Technical  pfote. 


A.  PERFORMING  ORG  REPORT  NUMBER 

Technical  Note  I97H-4  J 

8^. CON  TRACT  OR  GRANT  NUMBE  ROJ 


( I'iJ  77/-  / 7 V V \ /5[^ \ \9f,2H~-iH-('-yim 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 
Lincoln  Laboratory,  M.l.T.  1 
P.O.  Box  73 
Lexington,  MA  02173 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 
Air  Force  Systems  Command,  USAF 
Andrews  AFB 
Washington,  DC  20331 


14.  MONITORING  AGENCY  NAME  B ADDRESS  0/  different  front  ControUtmPlJice) 

Fleet ronic  Systems  Division  ( / J J nj-f r\ 

Han  scorn  AFB  ' i j ‘ 

Bedford,  MA  017.31  / I 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  8,  WORK  UNIT  NUMBERS 

S'  j Program  l lrmvui  ‘ - 657051 

( J/C  / Pm  Mini  PWJ  6491,  / 


^JA,  REPORT  DATE  _ 

111  16  May  **7h  ‘ 

-»T  NUMBER  OF  PAGES 
146 

IS.  SECURITY  CLASS,  (of  ihn  report/ 
Unclassified 

1 So.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  (An  Report ) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abtlract  entered  in  Bloch  20,  tf  different  from  Report/ 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  reverie  tide  if  neceeeary  and  identify  by  block  number ) 


flywheel 

rotors 


hardware 

bearings 


attitude  control 
high  efficiency 


high-speed 


70.  ABSTRACT  (Continue  on  reverie  etde  if  neceteery  and  identify  by  block  number/ 


The  following  paper  studies  the  possibilities  of  using  high-speed  flywheels  with  magnetic 
bearings  in  communications  spacecraft.  Hardware  point  designs  have  been  completed  for  a 
very  low  drag  bearing  and  a very  efficient,  compatible  motor  generator.  Using  the  results 
of  these  designs,  a tradeoff  study  compares  these  flywheels  with  combinations  of  NICd  batteries, 
N1H2  batteries,  and  conventional  low-speed  momentum  wheels. 


DD  , 1473  EDITION  OF  I NOV  *5  IS  OBSOLETE 


UNCLASS!  FI  p:D 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  Chen  Data  Entered ) 


<3o7  65~0 


13 


